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ABSTRACT 


Temperature  data,  collected  over  the  past  40  years,  was 
compiled  and  averaged  in  the  first  known  study  of  mean  ther- 
mal conditions  throughout  the  waters  of  Monterey  Bay.   The 
following  results  were  obtained: 

(1)   The  distribution  of  mean  sea  surface  temperature 
in  the  Bay  was  obtained  by  calculating  mean  monthly  values 
at  selected  grid  points  and  drawing  isotherms , 

(2)  Progressive  warming  of  the  upper  100m  was  observed 
to  occur  in  a  nearly  linear  fashion  with  time,  resulting  in 
a  maximum  temperature  increase  of  1.8°C  over  the  29  year 
interval  from  1931  to  1960.   This  implies  that  the  in- 
tensity of  upweiiing  is  also  decreasing  with  time, 

(3)  Geographic  variations  in  the  rates  of  upweiiing 
and  downwelling,  causing  a  relatively  "warm  tongue"  along 
the  Canyon  axis,  appeared  to  be  a  major  element  influencing 
the  dynamics  of  the  Bay,  and 

(4)  Mean  currents  were  inferred  from  the  topographies 
of  mean  monthly  isothermal  surfaces.   Limited  comparisons  of 
the  inferred  flow  to  past  measured  currents  were  very  en- 
couraging and  tended  to  support  the  feasibility  of  such  an 
inference  in  the  shallow  waters  and  irregular  topography  of 
Monterey  Bay. 
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I.   INTRODUCTION 

Monterey  Bay  is  a  large,  semi-elliptical  indentation 
in  the  California  coastline  separated  into  nearly  equal 
halves  by  a  deep  submarine  canyon  'Figure  1) . *   The  outer 
extremity  of  the  Bay  may  be  defined  by  a  line  joining  Point 
Santa  Cruz  to  the  north  and  Point  Pinos  to  the  south. 

A.    THERMAL  CONDITIONS  IN  MONTEREY  BAY 

Several  studies  have  been  conducted  of  thermal  condi- 
tions at  specific  locations  in  Monterey  Bay.   The  most  ex- 
tensive were  those  by  Skogsberg  [1936]  and  Bolin  [1964]  , 
both  spanning  five  year  periods.   Skogsberg  studied  the 
southern  portion  of  the  Bay  from  1929-1933  and  Bolin  re- 
searched conditions  at  a  station  in  the  Monterey  Submarine 
Canyon  (see  Appendix  A) .   Although  these  studies  well  de- 
fined the  thermal  conditions  for  a  portion  of  the  Bay,  ap- 
parently no  one  had  investigated  the  Bay  as  a  whole. 

As  discussed  in  Appendix  A,  the  author  compiled  a 
considerable  amount  of  temperature  (and  some  salinity)  data 
covering  the  period  of  time  from  Skogsberg 's  study  until 
the  present.   As  outlined  by  the  author's  objectives  (Para- 
graph C) ,  this  data  was  utilized  to  study  mean  monthly  ther 
mal  conditions  throughout  the  Bay  and  to  investigate  the 
possible  relationship  between  the  thermal  conditions  and 
indicated  mean  currents. 


i 

All  figures  will  be  found  in  Appendix  C 
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B.    CURRENTS  IN  MONTEREY  BAY 

1 .  Importance 

The  distribution  of  current  in  Monterey  Bay  is 
virtually  unknown  but  there  is  little  doubt  of  its  great 
importance  in  determining  the  ecological  future  of  the  area. 
Pollutants,  transported  by  the  currents,  are  scattered 
freely  along  the  beaches  in  the  form  of  sewage  and  general 
rubbish.   Less  evident  is  the  adverse  effect  on  marine  life 
created  by  the  discharge  of  industrial  pollutants  into  the 
Bay.   Currents  not  only  distribute  "harmful"  pollutants  but 
also  serve  to  maintain  the  chemical  and  thermal  balance 
necessary  to  support  marine  life. 

2 .  Why  Little  is  Known  Concerning  Distribution  of  Bay 
Currents 

The  well  established  currents  of  the  open  oceans 
are  relatively  strong  and  unchanging  with  time  (e.g.,  the 
Gulf  Stream  and  the  offshore  California  Current)  and  hence, 
their  general  movements  are  fairly  well  understood.   Monterey 
Bay,  being  enclosed  on  three  sides,  presents  a  different 
situation.   The  Bay  currents  are  weak  and  variable  and  their 
mean  flow  is  almost  completely  disguised  by  the  effects  of 
local  winds,  tides,  topography,  and  seasonal  changes  in  the 
offshore  currents . 

Past  current  studies  in  the  Bay  have  been  of  a 
limited  nature,  both  in  time  and  in  geographical  coverage. 
They  have  contributed  towards  an  understanding  of  "short 
term"  currents  (such  as  those  resulting  from  the  tides)  but 
too  few  observations  have  been  made  to  be  of  much  use  in 
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describing  the  mean  flow.   In  time,  and  with  many  addi- 
tional current  observations,  a  stud/  of  the  combined  ef- 
fects could  prove  to  be  of  great  value. 

3.   The  Feasibility  of  Approximating  Mean  Currents  by 
the  Geostrophic  Approach 

Indirect  approximations  of  current  using  the  geo- 
strophic theory  may  be  divided  into  two  categories:   (1) 
the  classical  use  of  the  geostrophic  method,  and  (2)  ap- 
proximations derived  from  the  classical  method. 

a.   The  Classical  Geostrophic  Method 

Basically,  flow  obtained  by  classical  applica- 
tion of  the  geostrophic  method  is  that  which  results  from 
the  balancing  of  the  Coriolis  force  and  the  pressure  force 
between  two  points  on  a  surface  of  variable  elevation, 
(Neumann  and  Pierson  [ly66]   discuss  the  theory  in  their 
book  entitled,  Principles  of  Physical  Oceanography) .   Com- 
puted elevation  of  the  sea  surface  at  a  given  location  in 
the  ocean  is  represented  by  a  term  called  the  "dynamic 
height,"  a  relative  parameter  based  on  the  density  of  sea 
water  at  the  location  and  referred  to  an  assumed  level-of- 
no-motion.   A  high  value  of  density  (usually  colder  water)' 
yields  a  lower  dynamic  height  and  conversely.   By  geo- 
strophic theory,  the  circulation  in  the  Northern  Hemisphere 
is  counterclockwise  about  locations  of  lower  dynamic  height 
(depressed  surfaces)  and  clockwise  about  locations  of  higher 
dynamic  heights  (elevated  surfaces).   The  level-of -no-motion 
is  the  depth  of  water  (as  best  can  be  ascertained)  where 
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absolute  current  velocity  equals  zero  and  this  rarely  oc- 
curs at  depths  less  than  150m. 

Accurate  representation  of  the  mean  current 
based  on  the  classical  geostrophic  method  thus  requires  two 
prerequisites;  (1)  a  knowledge  of  mean  temperature  and 
salinity  (to  obtain  mean  density)  of  the  water  column  at 
several  well  spaced  locations  in  the  area  of  interest,  and 
(2)  a  level-of -no-motion  everywhere  shallower  than  the 
depth  of  the  bottom  (otherwise,  there  is  no  reference  ex- 
cept the  bottom  itself  from  which  to  compute  true  velocity) 
Finally,  of  course,  wind  driven  and  tidal  currents  are 
neglected  in  computing  means. 

b.   Approximations  Derived  from  the  Classical 
Method 

Several  approximations  to  the  classical  ap- 
proach have  been  used  in  relatively  deep  waters  with  con- 
siderable success.   Barker  and  Denham  [1970]  used  the 
classical  method  to  plot  geopotential  topographies  (sur- 
faces of  equal  dynamic  height)  for  various  levels -of -no- 
motion  off  the  northeast  coast  of  New  Zealand.   The  close 
correlation  between  the  shapes  of  the  topographies  sug- 
gested that  broad  current  patterns  could  be  approximated 
without  using  an  accurate  level-of -no-motion .   Duncan  and 
Nell  [1969]  demonstrated  that  isopycnal  (equal  density) 
surfaces  were  closely  related  to  current  patterns  in  their 
study  off  the  Cape  Coast.   Drift  cards  were  employed  to 
measure  general  current  trends  and  were  found  to  follow  the 
isopycnal  contours  as  predicted  by  geostrophic  theory  for 
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the  Southern  Hemisphere.   Leipper  [1970]  demonstrated  that 
temperature  data  alone  could  be  used  to  approximate  cur- 
rent in  the  Gulf  of  Mexico.   Monthly  topographies  of  the 
22°C  isothermal  surface  were  constructed  by  him  from  data 
collected  on  various  cruises  during  the  period  1965-1966. 
The  current  was  found  to  closely  follow  the  contours,  moving 
in  a  direction  placing  colder  waters  to  the  left  of  the 
flow.   These  results  were  verified  by  geostrophic  flows  com- 
puted at  various  locations  in  the  Gulf.   This  approach  is 
based  on  the  assumption  that  temperature  changes  are  di- 
rectly proportional  to  density  changes  throughout  the  water 
column,  and  thus  isothermal  surfaces  may  be  considered  in 
place  of  isopycnal  surfaces. 

Application  of  any  of  the  above  approximations 
to  the  shallow  waters  and  irregular  topography  of  Monterey 
Bay  would  be  a  considerable  extension  of  usual  practice. 
However,  the  theory  of  geostrophic  flow  does  not  depend  on 
the  depth  of  water,  therefore,  geostrophic  methods  might  be 
useful  in  indicating  broad  current  patterns  in  the  shallow 
waters  of  the  Bay.   Considerable  temperature  data  was  com- 
piled for  Monterey  Bay  in  this  study  (see  Appendix  A) .   The 
author  used  this  data  to  investigate  the  feasibility  of 
relating  mean  currents  to  mean  isothermal  surfaces  in  a 
manner  similar  to  that  of  Leipper. 

C.    OBJECTIVES 

The  objectives  of  this  study  are  to: 
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1.  Compile  temperature  and  salinity  data  from  previous 
oceanographic  studies  conducted  in  Monterey  Bay, 

2.  Conduct  a  study  of  the  distribution  of  mean  monthly 
temperature  in  Monterey  Bay  and  construct  charts  of  the  mean 
monthly  distribution  of  sea  surface  temperature, 

3.  Investigate  the  feasibility  of  inferring  mean  cur- 
rents from  the  topographies  of  mean  isothermal  surfaces,  and 

4.  Prepare  a  computer  program  for  the  continuing 
storage  and  retrieval  of  oceanographic  data  collected  in 
Monterey  Bay. 

D.    ORGANIZING  AND  USING  THE  DATA 

The  IBM  360  computer  at  the  Postgraduate  School  was 
used  to  compile  the  temperature  data  and  to  compute  monthly 
averages  (as  discussed  in  Appendix  B) .   A  grid  system, 
dividing  the  Bay  into  19  blocks,  4mi  by  4mi  in  dimensions, 
was  used  to  segregate  the  temperature  data  (Figure  1) .   The 
author's  initial  grid  system  covered  a  larger  area  (Figure  2) 
containing  39  blocks,  but  only  the  numbered  blocks  held 
significant  data.   The  mean  temperature  of  a  given  block, 
based  on  the  averaging  of  all  temperatures  compiled  for  that 
block,  was  assumed  to  be  effectively  located  at  its  center. 

The  temperature  data  was  compiled  over  the  40  year 
interval  from  1929  through  1968.   Temperatures  collected  by 
all  reliable  methods  were  combined  in  the  averaging  process 
(although  the  computer  program  did  allow  for  data  to  be 
retrieved  by  type,  i.e.,  Nansen,  XBT,  BT) . 
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The  risks  inherent  in  the  "lumping"  together  of  all 
types  of  temperature  data  spanning  a  40  year  interval  and 
further  reducing  data  from  4m i  by  4r.ii  blocks  to  represent 
points  in  their  centers  were  fully  realized.   This  lumping 
was  essential  due  to  the  erratic  diitribution  of  data  within 
the  Bay  and  also  because  no  reasonably  short  time  span  held 
enough  data  to  support  this  study.   The  errors  introduced 
by  the  above  procedures  are  discussed  in  Chapter  II. 
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II.   VALIDITY  OF  THE  COMPUTED  MEAN  MONTHLY  TEMPERATURES 
(AND  DENSITIES) 

A  discussion  of  the  questions  raised  in  the  Introduc- 
tion involving  the  method  of  computing  mean  monthly  temper- 
atures is  considered  necessary  before  an  appraisal  of  the 
results  (Chapters  III  and  IV)  can  be  made. 

A.    TIME  RANGE  OF  DATA 

Data  from  three  separate  time  intervals  within  the 
previously  mentioned  40  year  time  span  was  averaged  and 
compared  to  determine  the  degree  of  consistency  of  thermal 
conditions  relative  to  the  progression  of  time.   Three  Bay 
studies  (as  referenced  in  Appendix  A)  were  used:   (1) 

CalCOFI  surveys  (1954- 1967)  . 2   The  time  intervals  between 
midpoints  (mean  years)  of  the  studies  were  22  years  from 
Skogsberg  to  Bolin  and  7  years  from  Bolin  to  CalCOFI,  with 
a  total  time  interval  of  29  years  (1931  to  1960)  . 

The  annual  thermal  cycles  of  the  three  periods  (Fig- 
ures 4-6)  were  generally  similar  in  shape  with  upwelling 
noticeable  from  the  surface  to  500m  (as  characterized  by 
the  humps  in  the  isotherms  during  the  months  of  April 
through  June) .   A  consistent  feature  over  the  29  year  in- 
terval was  the  apparent  weakening  in  intensity  of  upwelling 


2 

The  oceanographic  stations  used  in  these  studies  are 
shown  in  Figure  3.   Skogsberg 's  station  is  denoted  by  "C", 
Bolin' s  by  "B"  and  CalCOFI  by  a  square  numbered  "3".   It  is 
noted  that  the  three  stations  are  contained  within  a  radius 
of  two  mi . 
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with  time.   The  9°C  isotherms  from  all  three  studies  are 
plotted  in  Figure  4  to  demonstrate  the  effect.   It  was  ob- 
served that  the  depths  of  the  9°C  isotherms  were  relatively 
close  together  during  the  winter  months.   During  the  re- 
mainder of  the  year  (especially  noticeable  in  May)  the  9°C 
isotherms  were  observed  to  deviate  from  one  another.   The 
deviation  followed  a  consistent  pattern  with  the  9°C  iso- 
therm of  the  earliest  study  (Skogsberg  in  1931)  at  the 
shallowest  depth  and  the  9°C  isotherms  of  the  later  two 
studies  at  greater  depths  as  time  progressed.   This  "flat- 
tening out''  of  the  9°C  isotherm  with  time  could  indicate  a 
weakening  in  intensity  of  upwelling. 

A  general  warming  with  time  was  also  observed  to  occur 
throughout  the  waters  of  the  Bay.   Figure  6A  depicts  the 
temperature  increase  (AT)  observed  to  occur  over  the  period 
from  1931  (Skogsberg's  mean)  to  1960  (CalCOFI's  mean)  at 
selected  depths  throughout  the  year.   The  curves  represent 
the  variation  of  AT  for  each  month  and  for  depths  from  the 
surface  to  500m.   Figure  6A  shows  that  the  maximum  tempera- 
ture increase  occurred  at  depths  from  the  surface  to  100m 
and  during  the  period  from  October  to  December.   At  depths 
from  200-500m,  AT  cycled  between  positive  and  negative 
values  throughout  the  year;  however,  AT  was  always  positive 
(temperature  increase)  in  the  months  of  May,  June,  September, 
and  October. 

The  consistently  positive  values  of  AT  in  May  and  June 
for  depths  from  the  surface  to  500m  might  be  another  indica- 
tion that  upwelling  was  decreasing  in  intensity  (with  time) 
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throughout  this  region,  since  the  warming  in  May  and  June 
would  also  indicate  a  flattening  out  of  the  isotherms. 
The  reasons  for  the  maximum  temperature  increase  during  the 
months  of  October  to  December  are  unknown.   It  could  be  re- 
lated to  changes  in  the  "Oceanic  Period,"  or  possibly,  to 
changes  in  the  transition  period  between  the  influence  of 
the  California  and  Davidson  current  systems. 

In  summary,  the  general  warming  of  the  water  column 
with  time  was  observed  to  be  consistent  in  the  upper  100m 
(by  the  similarity  in  the  shape  of  the  AT  curves  in  Figure 
6A  at  these  depths).   At  depths  below  200m,  the  consistency 
was  not  as  good  (as  observed  by  the  AT  curves  cycling  be- 
tween positive  and  negative  values);  however,  the  total 
range  of  AT  was  only  0.5°C  in  the  deeper  waters. 

The  maximum  temperature  increase  over  the  29  year 
interval  was  1.8°C  (the  50m  curve  in  November  from  Figure 
6A)  resulting  in  about  l/15th  of  a  degree  C  increase  per 
year.   This  small  annual  temperature  change,  along  with 
the  general  similarity  of  the  curves  (as  observed  in  Fig- 
ures 4-6)  was  considered  to  represent  consistent  behaviour 
with  time.   It  was  concluded  that  the  author's  method  of 
averaging  temperature  values  was  justified.   This  provided 
averages  that  were  representative  of  the  overall  period. 

B.    POSITIONING  OF  STATIONS 

The  computer  grid  system  was  based  on  the  assumption 
that  all  temperature  values  found  in  a  given  block  were  lo- 
cated at  its  center.   The  grid  system  was  initially  positioned 
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so  as  to  minimize  the  distance  of  known  oceanographic  sta- 
tions from  grid  centers.   The  maximum  deviation  from  such 
centers  was  1.6  mi.   Data  was  assumed  to  be  scattered 
randomly  around  the  centers  of  the  blocks  in  the  grid  sys- 
tem. 

C.    USE  OF  THREE  DIFFERENT  DATA  COLLECTION  SYSTEMS 

The  data  was  collected  by  Nansen,  BT ,  and  XBT  devices. 
Several  papers  have  been  written  concerning  system  errors 
peculiar  to  each.   The  magnitude  of  the  errors  in  the  BT ' s 
and  XBT ' s  are  apparently  undeterminable. 

To  investigate  the  significance  of  possible  system 
errors,  ths  standard  deviations  (a)  of  temperature  (and 
density)  values  were  included  in  the  computer  output. 
(X  ±  a)  represented  the  interval  containing  671  of  all 
values  averaged  for  a  given  mean  (X) . 

Standard  deviations  of  mean  monthly  sea  surface  temper- 
ature versus  time  (by  months)  were  plotted  for  selected 
blocks.   Figure  7  (Block  13)  represented  the  worst  condi- 
tion found.   The  larger  standard  deviation  in  July  was  not 
representative;  the  consistent  features  of  all  graphs  were: 
(1)  minimum  values  of  standard  deviation  for  the  month  of 
January,  and  (2)  the  relatively  small  variation  of  the 
values  throughout  the  year. 

Plots  of  standard  deviation  versus  depth  (Figure  8) 
were  also  consistent.   Values  of  standard  deviation  halved 
within  30 -40m  of  the  surface,  reached  values  of  0.5°C  anywhere 
from  50-275m  and  then  linearly  decreased  to  minimum  values  of 
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0.1-0.2°C  in  deeper  waters.   The  decrease  of  standard  devia- 
tion with  increasing  depth  of  water  and  its  practically- 
invariable  value  (0.1-0.2°C)  below  800m,  where  the  data  was 
exclusively  from  Nansen  casts,  indicated  that  seasonal 
changes  were  mainly  responsible  for  variations  in  mean 
temperatures.   System  errors  were  thus  neglected. 

D.    GEOGRAPHIC  DISTRIBUTION  OF  DATA 

The  3000  oceanographic  stations  used  in  this  study 
(25,100  temperatures)  were  erratically  distributed  through- 
out the  blocks  in  the  computer  grid  system  (Figure  2) . 
Blocks  1-6,  13,  20,  and  21S  contained  the  major  portion  of 
the  data  (10-85  surface  temperatures  per  block  per  month) 
and  thus  were  given  major  consideration  in  the  analysis. 
The  remainder  of  the  blocks  contained  random  amounts  of 
data,  except  for  Blocks  17  and  19  which  generally  contained 
3-6  values  per  month  each. 

Clearly,  the  blocks  containing  relatively  few  tempera- 
tures (less  than  10  per  month)  would  yield  less  reliable 
averages  than  blocks  with  many  values.   Choosing  a  minimum 
number  of  temperature  values  to  achieve  relatively  accurate 
mean  monthly  temperatures  was  based  on  a  study  of  continuity 
between  weak  blocks  (with  less  than  10  temperatures  per 
month  at  a  given  depth)  and  strong  blocks  (more  than  10  tem- 
peratures) . 

Graphs  of  annual  mean  temperature  cycles,  mean  monthly 
temperature  sections,  and  standard  deviation  versus  time  and 
depth  were  drawn  for  weak  Blocks  8,  17,  and  19  and  for  strong 
Blocks  1-6  and  13.   The  annual  mean  temperature  cycles  of  the 
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weak  blocks  compared  closely  to  those  of  the  strong  blocks. 
This  comparison  was  well  represented  by  Block  19  (Figure  9) 
and  Block  3  (Figure  10)  . 

Further  continuity  in  data  between  weak  and  strong 
blocks  was  demonstrated  by  Figure  11,  showing  the  thermal 
structures  of  two  vertical  sections  (see  Figure  12  for  lo- 
cation of  sections) ,  one  oriented  north  to  south  connecting 
Blocks  1,  2,  17,  and  5  (S-l)  and  the  other  oriented  along 
the  Monterey  Submarine  Canyon  axis  (S-2)  connecting  Blocks 
13,  3,  17,  and  19.   The  isotherms  were  smooth  and  continuous 
in  the  upper  and  lower  waters  for  both  sections  throughout 
the  year.   In  depths  from  50-100m,  the  isotherms  were  smooth 
and  continuous  along  section  S-2.   Along  section  S-l  the  iso- 
therms were  occasionally  depressed  at  Block  17.   The  depres- 
sion of  isotherms  at  Block  17  in  section  S-l  and  lack  of 
depression  for  the  same  block  in  section  S-2  will  be  shown 
to  be  consistent  in  Chapter  III,  paragraph  B. 

Finally,  the  standard  deviations  of  mean  monthly  tem- 
peratures versus  time  and  depth  were  comparable  for  all 
blocks,  weak  or  strong  (Figure  7  and  8  were  representative). 

The  continuity  exhibited  in  the  vertical  sections  be- 
tween adjacent  weak  and  strong  blocks  (and  within  the  blocks 
as  time  progressed)  supported  the  limited  use  made  of  the 
weak  blocks  for  constructing  mean  monthly  isothermal  sur- 
faces in  Chapters  III  and  IV.   As  the  isothermal  surfaces 
were  prepared,  the  weak  blocks  demonstrated  further  consis- 
tency of  data,  especially  Blocks  8,  17,  and  19.   Examples  of 
this  are  noted  in  Chapter  IV. 
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III.   THERMAL  CONDITIONS 

A.    THE  ANNUAL  THERMAL  CYCLE  IN  MONTEREY  BAY 

The  "3-phase"  annual  thermal  cycle  occurring  at  specific 
locations  in  Monterey  Bay  has  been  described  by  Skogsberg 
[1936],  Bolin  [1964]  and  others.   The  following  discussion 
considers  the  variation  in  mean  thermal  conditions  relative 
to  various  geographic  locations  throughout  the  Bay,  a  sub- 
ject apparently  not  previously  dealt  with. 

Figure  12  depicts  three  geographic  regions  in  the  Bay 
that  the  author  will  refer  to  in  this  Chapter.   The  choice 
of  the  50  fathom  curve  as  a  dividing  line  was  arbitrary,  but 
gave  the  desired  sectioning. 

Tno      f  onirnn       Doni  <~>-r» 

Since  Bolin  described  the  thermal  conditions  for 
the  Canyon  region  a  detailed  account  will  not  be  repeated 
here;  however,  a  brief  description  of  the  "phases"  will  be 
necessary  for  later  comparisons. 

a.   The  "Upwelling  Period" 

Commonly  considered  to  span  from  mid-January  to 
September,  the  Upwelling  Period  is  characterized  by  the 
rising  of  the  deep,  colder  waters.   As  found  in  this  study 
(Figure  10) ,  the  cold  waters  began  their  upward  movement  in 
mid- January  from  depths  as  great  as  700m  (however,  the  move- 
ments were  minimal  below  300m) .   The  maximum  rate  of  rise 
occurred  between  February  and  March  and  the  peak  was  found 
in  April  or  May  when  the  11°C  water  nearly  reached  the 
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surface.   Although  upwelling  reached  its  peak  later  at 
greater  depths,  none  was  observed  after  July. 

The  descent  (downwelling)  of  the  cold  waters 
generally  commenced  in  June  and  was  apparent  until  late 
November  at  greater  depths. 

b.  The  "Oceanic  Period" 

This  period  is  commonly  considered  to  occur 
during  September  and  October  and  is  well  defined  by  the 
sharp  vertical  temperature  gradients  at  depths  from  the 
surface  to  about  100m. 

c.  The  "Davidson  Current  Period" 

In  November  the  south-moving  California  Cur- 
rent is  displaced  by  the  north-moving  Davidson  Current  along 
the  coast.   The  influx  of  warmer  waters  into  the  Bay  re- 
sults in  a  weakening  of  the  vertical  temperature  gradients 
and  a  resulting  well-mixed  layer  of  water  from  the  surface 
to  about  50m.   The  northerly  flow  continues  until  mid- 
January  when  the  California  Current  again  predominates  and 
the  "3-phase"  annual  cycle  commences  again. 

2 .   The  Regions  of  the  Shallows  (see  Figure  12) 

The  annual  thermal  cycle  in  the  Shallows  was  simi- 
lar to  that  in  the  Canyon  Region  with  two  noteworthy  ex- 
ceptions : 

First,  the  very  shallow  regions  at  the  extreme 
northern  and  southern  ends  of  the  Bay  possessed  unique 
annual  thermal  cycles  at  depths  from  the  surface  to  10m. 
Referring  to  Figure  13,  the  annual  cycle  of  the  surface 
temperature  in  the  Canyon  (Block  3)  showed  a  clear  decrease 
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in  temperature  from  January  until  early  May  and  then  a 
steady  increase  until  late  September  (this  was  typical  for 
all  the  deep  water  areas) .   The  annual  cycles  of  surface 
temperature  for  the  very  shallow  waters  (Blocks  1  and  5) 
showed  no  similar  decrease  in  late  winter  but  rather  a 
steady  increase  from  mid-January  until  late  September 
(Figures  13  and  14).   This  steady  temperature  increase, 
unique  for  the  shallow  water  regions  (less  than  75m) ,  was 
probably  due  to  greater  heating  of  the  water  column  from 
bottom  reflection  and  the  insolation  of  nearby  land  masses. 
This  greater  heating  of  the  shallow  areas  is  clearly  shown 
in  the  distribution  of  mean  monthly  sea  surface  temperatures 
as  will  be  discussed  in  paragraph  E  (Figures  19-30).   The 
relatively  deep  regions  of  the  Shallows  areas  (Blocks  2  and 
4)  exhibited  a  very  slight  amount  of  cooling  from  January 
until  May,  but  not  nearly  as  much  as  did  the  can/on  Region 
(see  Figures  15  and  16).   The  second  exception  is  discussed 
below. 

B.    THE  VARIATIONS  IN  INTENSITY  OF  UPWELLING 

In  studying  the  various  graphs  constructed  for  the 
analysis  of  thermal  conditions,  the  author  came  upon  a 
feature  which  could  be  an  important  element  in  the  dynamics 
of  the  Bay  in  general.   Figure  17  is  introduced  to  show  the 
months  of  the  maximum  rates  of  upwelling  and  downwelling 
during  the  year.   Inspection  of  the  annual  cycle  in  Block  3 
shows  maximum  rates  of  upwelling  from  February  to  March,  and 
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maximum  rates  of  downwelling  (or  heating)  from  July  to 
August  (note  the  slopes  of  the  10°C  and  11°C  isotherms). 

Vertical  temperature  section  S-l  (see  Figure  12)  shows 
a  strange  dipping  of  the  isotherms  at  Block  17  for  all  of 
the  above  mentioned  months  (Figure  18  is  representative) . 
Block  17  was  a  weak  block  (less  than  10  temperatures  per 
month)  and  thus  the  dipping  was  suspected.   To  attempt  to 
verify  the  dipping  of  isotherms  in  the  Canyon,  two  addi- 
tional sections  were  constructed  (see  Figure  12  for  loca- 
tions) ,  both  running  across  the  Canyon  axis  (as  did  section 
S-l).   Section  S-3  connected  Blocks  2,  3,  and  4,  and  sec- 
tion S-4  connected  Blocks  2,  3,  and  6  (emphasis  is  placed 
on  the  fact  that  sections  S-3  and  S-4  contained  only  strong 
blocks  with  greater  than  20  temperatures  per  month  each) . 3 
The  results  are  clearly  illustrated  in  Figure  J  8  which  shows 
marked  dipping  of  the  isotherms  across  the  Canyon  axis  in 
all  sections.   This  results  in  a  relatively  warm  tongue 
along  the  Canyon  axis  for  the  months  of  maximum  rates  of  up- 
welling  and  downwelling. 

The  dipping  of  the  isotherms  across  the  Canyon  axis 
could  have  been  caused  by  the  up- slope  entry  of  the  cold  up- 
welling  waters  into  the  regions  of  the  Shallows.   During 
February  and  March,  when  upwelling  rates  are  highest,  the 
cold  waters  are  ascending  from  depths  down  to  300m.   As  the 


3 

Table  V  lists  the  number  of  temperatures  compiled  fo: 
the  various  blocks  of  the  computer  grid  system. 
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cold,  dense  water  spreads  to  the  Canyon  walls,  its  hori- 
zontal movement  is  interrupted,  and  thus  it  is  accelerated 
up  the  walls  and  into  the  shallower    adjacent  regions. 
Conversely,  in  July  and  August,  when  downwelling  rates  are 
highest,  the  cold  waters  are  descending  from  the  near- 
surface  regions.   Their  rates  of  descent  are  accelerated 
into  the  Canyon  since  the  deepening  topography  offers  no 
resistance  to  their  sinking.   The  transition  period  (from 
April  to  June)  between  upwelling  and  downwelling  showed  no 
dipping  of  isotherms  across  the  Canyon. 

No  dipping  of  isotherms  was  observed  along  the  Canyon 
axis.   This  would  be  as  expected  since  the  depth  of  water 
along  the  axis  (section  S-2)  was  everywhere  deeper  than  300m 
thus  affording  no  bottom  influence  unon  the  cold  uowelling 
waters . 

The  annual  mean  temperature  cycles  of  the  relatively 
deep  regions  (Figure  16)  and  the  very  shallow  regions  (Fig- 
ure 14)  of  the  Shallows  further  demonstrated  the  up-slope 
surge  idea.   During  the  months  of  highest  up/downwelling 
rates,  Blocks  2  and  4  (adjacent  to  the  Canyon)  exhibited 
somewhat  stronger  upwelling  (greater  slopes  of  the  10°C 
and  11°C  isotherms)  than  did  Blocks  1  and  5  which  were  at 
the  extreme  northern  and  southern  ends  of  the  Bay.   It  thus 
appeared  that  the  cold  up-surging  waters  might  have  lost 
some  of  their  momentum  by  the  time  they  reached  the  very 
shallow  regions  of  Blocks  1  and  5. 
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C.    COMPARISON  OF  THE  ANNUAL  THERMAL  CYCLE  IN  THE  BAY 
TO  THAT  OF  AN  OFFSHORE  LOCATION 

Block  20,  located  22  mi  west-southwest  of  the  Bay 

(Figure  2)  was  compared  to  Block  3  in  the  Canyon.   Figure 

17  demonstrated  two  broad  differences  between  the  thermal 

conditions  of  the  two  locations. 

1 .  Mixing  in  the  Upper  Layers 

The  offshore  station  exhibited  generally  warmer 
surface  waters  throughout  the  year  and  better  mixing  of  the 
upper  waters  (0-50m)  from  January  to  April.   At  onset  of 
the  "Oceanic  Period"  in  September,  the  offshore  and  Bay 
stations  both  exhibited  maximum  vertical  temperature  gradi- 
ents . 

At  the  close  of  the  "Oceanic  Period"  in  November, 
the  Bay  station  exhibited  somewhat  better  mixing  than  did 
the  offshore  station,  maintaining  its  more  uniform  surface 
waters  until  January.   The  weakening  of  the  vertical  tem- 
perature gradients  in  the  Bay  (as  described  by  Bolin)  was 
assumed  to  be  the  result  of  the  influx  of  warm  waters 
brought  into  the  Bay  by  the  north-flowing  Davidson  Current 
(during  the  period  of  November  through  January) .   The  more 
uniform  mixing  of  the  Bay  waters  might  have  been  due  to  a 
greater  influence  by  the  Davidson  Current  nearshore. 

2 .  Upwelling 

Upwelling  is  considered  to  result  from  the  northerly 
winds  deflecting  the  surface  waters  to  the  right  and  away 
from  the  coast,  allowing  the  deeper  cold  waters  to  rise.   The 
intensity  of  upwelling  should  then  decrease  as  the  distance 
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to  the  coast  increases.   Figure  17  demonstrated  that  this 
was  indeed  the  case,  Block  20  showing  less  intense  upwelling 
than  Block  3. 

D.  PROGRESSIVE  WARMING 

Chapter  II,  paragraph  A,  contains  a  discussion  of  the 
progressive  warming  found  for  the  entire  water  column  in 
the  Bay  over  the  29  year  interval  from  1931  to  1960.   It 
should  be  noted  that  Bolin  [1964]  found  an  opposite  indica- 
tion of  progressive  cooling  from  1951-1955.   He  investigated 
the  "previous  and  subsequent"  sea  surface  temperatures  taken 
at  selected  locations  along  the  west  coast  of  North  America 
and  found  that  the  years  of  his  study  were  all  relatively 
cool  and  that  1955  had  the  minimum  temperature  in  a  "long- 
term  trend."   In  Chapter  II,  evidence  was  presented  indica- 
ting continual  progressive  warming  over  the  29  year  interval 
noted  above.   It  thus  appears  that  che  cooling  observed  by 
Bolin  was  a  short-term  cycle. 

E.  THE  DISTRIBUTION  OF  MEAN  MONTHLY  SEA  SURFACE  TEMPERATURE 
IN  MONTEREY  BAY 

Figures  19-30  depict  the  distribution  of  mean  monthly 
sea  surface  temperature  in  Monterey  Bay. 

1 .   Description  of  Figures 

The  dots  on  each  figure  (representing  the  centers 
of  their  respective  blocks)  served  as  the  points  from  which 
isotherms  were  drawn.   The  values  of  mean  monthly  tempera- 
ture are  directly  below  the  points.   Table  I  lists  the  total 
number  of  temperature  measurements  per  block  per  month  (at 

29 


TABLE  I .   The  Range  in  the  Total  Number  of  Temperatures 
Averaged  to  Compute  the  Means  for  Blocks  in 
the  Bay 


BLOCK  DEPTH   RANGE  OF  TOTAL  TEMPERATURES  PER  MONTH  AT  0  5  50m 
NO.    (m)    J    F    M   A   M   J    J   A    S    0    N    D 
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Note  1:   Blocks  9-11,  14,  and  18  seldom  contained  data. 
Note  2:   Block  20  (see  Figure  3)  is  located  22  mi  offshore 
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the  surface  and  50m)  that  were  used  to  compute  the  mean 
temperatures.   Mean  temperatures  are  not  listed  (or  used) 
in  the  figures  if  less  than  four  values  were  available 
for  averaging.   Dashed  lines  represent  probable  extensions 
of  the  isotherms. 

The  circle  in  the  upper  right  hand  corner  of  each 
figure  displays  the  local  mean  wind  condition  (taken  from 
the  U.  S.  Naval  Oceanographic  Office  "Pilot  Charts," 
series  H.O.  1401)  and  the  mean  offshore  currents  (taken 
from  the  U.  S.  Naval  Oceanographic  Office  Publication  H.O. 
570) .   The  local  mean  winds  are  shown  as  a  series  of  ar- 
rows (fine  lines)  that  fly  with  the  winds.   The  lengths  of 
the  arrows  represent  the  percentages  of  the  month  that  winds; 
were  from  the  directions  shown  (the  percentage  shown  for 
the  longest  arrow  serves  as  a  reference) .   The  number  of 
tails  on  a  given  arrow  defines  the  wind  force  by  the 
Beaufort  Scale,  as  shown  below: 

Number  of  Tails       Mean  Wind  Force  (kt) 
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The  mean  offshore  current  is  displayed  by  an  arrow  (bold 
line)  extending  the  diameter  of  the  circle.   Mean  current 
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direction  is  with  the  arrow  and  its  speed  (miles  per  day) 
is  displayed  inside  the  inner  circle. 

The  upper  left  hand  corner  of  each  figure  shows  the 
mean  sea  surface  temperature  at  a  location  15  mi  offshore. 
This  data  was  taken  from  Naval  Oceanographic  Office  Special 
Publication,  SP-123,  and  represented  mean  monthly  values 
based  on  ship's  injection  temperatures  collected  over  the 
past  100  years . 

2 .   Discussion 

a.   The  Annual  Cycle  of  Mean  Sea  Surface  Temperature 
The  following  two  factors  appear  to  influence 
the  distribution  of  mean  sea  surface  temperature. 

(1)  Upwelling 

The  coldest  waters  in  the  Bay  occur  during 
the  period  from  February  to  May  (with  April  being  the  coldest 
month  of  the  year)  coinciding  with  the  upwelling  period. 
Figures  20-22  show  the  intrusion  of  the  colder  waters  from 
seaward.   It  appears  that  the  main  intrusion  broadly  oc- 
curs along  the  axis  of  the  Submarine  Canyon  (note  the  11°C 
and  12°C  isotherms  in  Figure  22) . 

(2)  The  Variation  in  Mean  Air  Temperature 
According  to  data  contained  in  the  U.  S. 

Naval  Oceanographic  Office  Pilot  Chart  series,  the  mean  air 
temperature  is  at  a  minimum  in  January,  rises  to  an  annual 
maximum  in  September  to  October,  and  then  falls  to  a  minimum 
again  in  January. 

The  relation  of  mean  sea  surface  to  mean 
air  temperatures  is  demonstrated  by  following  the  monthly 
position  of  the  13°C  isotherm.   In  January  and  February 
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absence  implies  generally  colder  waters.   In  March  it  was 
seen  to  intrude  into  the  Bay  near  the  Moss  Landing  area. 
In  April,  the  13°C  isotherm  appears  to  be  pushed  aside  by 
the  intrusion  of  the  colder  waters  (Figure  22) .   During  the 
period  from  May  to  October,  the  13°C  isotherm  was  observed 
to  move  steadily  seaward  as  the  waters  in  the  Bay  continued 
to  warm.   Although  the  positions  of  the  13°C  isotherms  are 
nearly  identical  in  September  and  October,  the  overall 
thermal  distribution  shows  September  to  be  the  warmest  month 
of  the  year  (Figure  27).   From  October  until  January,  the 
13°C  isotherm  was  observed  to  move  steadily  towards  shore 
indicating  the  broad  cooling  of  the  surface  waters. 

b.   Horizontal  Sea  Surface  Temperature  Gradients 
Horizontal  temperature  gradients  are  weakest 
during  the  period  from  November  through  January  with  nearly 
uniform  temperature  throughout  the  Bay  (less  than  0.5°C 
change)  occurring  in  January.   The  uniform  temperatures 
characterizing  this  period  are  probably  due  to  the  influx 
into  the  Bay  of  warmer  waters  from  the  north-flowing 
Davidson  Current  (prevailing  during  the  above  period) .   In 
February,  the  temperature  gradients  are  stronger  and  con- 
tinue to  increase  in  strength  until  reaching  their  annual 
maximum  in  July.   The  mean  offshore  sea  surface  temperature 
of  13.9°C  in  July  indicates  that  strong  gradients  are  also 
set  up  offshore  as  well  as  in  the  Bay.   Steady  weakening  of 
temperature  gradients  occur  from  August  until  January,  com- 
pleting the  annual  cycle. 
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It  will  be  recalled  that  maximum  vertical  tem- 
perature gradients  occur  in  the  surface  layers  during  the 
Oceanic  Period  from  September  to  October  (see  Figure  10). 
Bolin  related  the  occurrence  of  maximum  vertical  gradients 
to  the  relatively  windless  (and  also  the  minimum  velocity 
of  the  California  Current)  conditions  prevalent  during  this 
period.   The  occurrence  of  the  maximum  horizontal  sea  sur- 
face temperature  gradients  in  July  (2-3  months  earlier  than 
the  maximum  vertical  gradients)  might  be  related  to  the 
greater  velocity  of  the  south- flowing  California  Current. 
The  cooler  waters  carried  into  the  Bay  by  this  current  move 
at  twice  the  velocity  in  July  as  compared  with  September 
to  October  (3.6  versus  1.7  mi/day).   Hence,  the  outer  Bay 
regions  are  colder  whi]e  the  inner  or  nearshore  regions 
continue  to  warm  from  atmospheric  heating. 

The  mean  sea  surface  temperatures  in  the  north- 
ern and  southern  areas  of  the  Bay  are  consistently  warmer 
than  those  in  the  center  area  for  all  months  except 
November  through  January  (when  generally  uniform  tempera- 
tures prevail) .   The  northern  area  exhibited  consistently 
higher  temperatures  than  did  the  southern  area  throughout 
this  period.   The  warm  spot  in  the  northern  region  is  pro- 
bably due  to  the  shallower  depths  of  water  there  (see 
Figure  12)  thus  permitting  greater  insolation  and  reflective 
heating  of  the  waters  as  earlier  mentioned.   Also,  the  shal- 
lower northern  region  is  somewhat  isolated  from  the  effects 
of  upwelling  and  thus  does  not  undergo  mixing  with  the 
deeper  and  colder  waters. 
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3 .   Summary 

Inspection  of  the  mean  isothermal  surfaces  por- 
trayed in  Figures  19-30  reveal  several  blank  spots  in  the 
data  field,  especially  in  the  northeastern  section  of  the 
Bay.   It  would  be  of  considerable  benefit  to  conduct  a 
further  search  for  data  within  these  blank  spots  (Blocks 
16,  18,  and  the  non-numbered  block  shown  crossed  out  in 
Figure  1)  thus  providing  a  more  complete  picture  of  the 
sea  surface  temperature  distribution. 
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IV.   USE  OF  ISOTHERMAL  SURFACES  AS  INDICATORS  OF  MEAN  CURRENT 

If  mean  flow  can  be  approximated  in  Monterey  Bay  using 
mean  isothermal  surfaces,  then  the  topographies  of  these 
surfaces  must  be  associated  with  all  dynamic  forces  acting 
to  establish  mean  flow.   Some  of  the  factors  involved  are 
as  follows : 

A.    POSSIBLE  GENERATING  FORCES 
1 .   Offshore  Currents 

The  offshore  currents  are  considered  by  some  to  be 
the  major  driving  forces  behind  the  generation  of  Bay  cur- 
rents. Some  theorists  believe  that  a  northerly  flow  adja- 
cent to  the  Bay  will  establish  a  clockwise  gyre  within  it, 
and  conversely  for  a  southerly  flow.  This  theory  was  pre- 
sented recently  by  Garcia  [1971]  using  a  mathematical  model 
of  the  Bay. 

a.  The  California  Current 

Transgressing  a  large  clockwise  circle  in  the 
Northern  Hemisphere,  the  California  Current  makes  its  way 
south  along  the  coast  year-round.   The  southern  flow  of  its 
outer  boundaries  is  felt  at  the  seaward  edge  of  the  Bay  from 
mid-February  until  early  November,  with  mean  velocities 
ranging  from  1.0  to  6.7  mi/day  (see  Chapter  III,  paragraph 
E.l). 

b.  The  Davidson  Current 

Prevailing  from  November  to  early  February,  the 
Davidson  Current  flows  north  adjacent  to  the  coast,  bringing 

relatively  warmer  waters  into  the  Bay  (see  Figure  31) . 
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2.  Winds 

Long  acting  and  relatively  steady  winds  will 
generate  related  current  patterns  wherever  and  whenever 
such  wind  systems  exist. 

3.  Tides 

Tidal  currents  are  active  in  Monterey  Bay  but  are 
not  considered  in  this  study  since  they  are  very  nearly 
(if  not  entirely)  averaged  out  in  taking  monthly  means. 

4.  Upwelling 

Upwelling  has  been  studied  as  thoroughly  as  any 
oceanographic  feature,  but  almost  all  studies  considered 
only  its  vertical  aspects.   As  discussed  in  Chapter  III, 
paragraph  B,  the  author  found  that  the  horizontal  dif- 
ferences in  upwelling  in  the  vicinity  of  the  Submarine  Can- 
yon were  apparently  a  contributing  factor  in  the  resulting 
shapes  of  the  mean  isothermal  surfaces  (possibly  a  major 
factor  during  the  months  of  maximum  rates  of  up/downwelling) . 
The  possible  influence  of  up/downwelling  upon  mean  indicated 
currents  is  discussed  in  paragraph  D  of  this  chapter. 

B.    VERIFYING  THE  CORRELATION  BETWEEN  ISOTHERMAL  AND  ISO- 
PYCNAL  SURFACES 

To  employ  an  approach  similar  to  that  used  by  Leipper 
[1970]  in  the  Gulf  of  Mexico  (discussed  in  the  Introduction), 
it  was  first  necessary  to  show  that  changes  in  mean  density 
were  directly  proportional  to  changes  in  mean  temperature 
throughout  the  water  column  in  the  Bay.   There  was  not 
enough  data  available  on  mean  density  to  allow  the  con- 
struction of  vertical  density  sections,  therefore,  annual 
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cycles  of  mean  density  and  temperature  (for  several  blocks) 
were  compared.   The  results  showed  that  mean  density  and 
temperature  were  closely  related  throughout  the  Bay  and  at 
all  depths  considered  except  in  the  region  above  10m  where 
the  correlation  was  not  quite  as  good.   The  below  listed 
comparisons  demonstrated  the  results: 

(1)  Figure  32  (mean  density  structure,  Blocks  1  and  5)  to 
Figure  14  (mean  temperature  structure,  Blocks  1  and  5), 

(2)  Figure  32  (mean  density  structure,  Block  20)  to 
Figure  17  (mean  temperature  structure,  Block  20), 

(3)  Figure  33  (mean  temperature  structure  with  mean  density 
structure  (dashed  lines)  superimposed,  Block  13). 

The  relatively  small  total  change  in  mean  density  through- 
out the  water  column  was  not  unexpected.   Previously  measured 
currents  in  the  Bay  were  generally  on  the  order  of  U.25  to 
0.5  kt ,  and  rarely  (if  ever)  exceeded  1  kt.   Therefore,  the 
small  changes  in  density  (assuming  isopycnal  surfaces  are 
related  to  indicated  mean  current)  would  imply  small  values 
of  mean  current.   This  appears  to  be  the  case. 

C.    PORTRAYAL  OF  THE  ISOTHERMAL  SURFACES 

The  foregoing  results  provide  evidence  demonstrating 
the  close  correlation  between  mean  isothermal  and  isopycnal 
surfaces.   Assuming  that  Leipper's  approximate  geostrophic 
approach  could  be  extended  from  deep  waters  to  the  waters  of 
Monterey  Bay,  the  mean  isothermal  surfaces  would  then  ef- 
fectively replace  the  mean  isopycnal  surfaces  and  indicated 
mean  currents  could  be  deduced  from  them.   On  this  basis, 
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topographies  of  the  mean  monthly  isothermal  surfaces  were 
assumed  to  represent  mean  flow  and  were  portrayed  in  Fig- 
ures 34-48. 

1 .  Selecting  the  Isothermal  Surface  for  Portrayal 
The  annual  variation  in  upwelling,  causing  the 

isotherms  to  rise  and  fall  through  a  considerable  depth  in- 
terval, prevents  the  use  of  any  one  isothermal  surface  for 
the  whole  year.   Instead,  isotherms  were  chosen  for  a  given 
month  by  inspection  of  selected  vertical  temperature  sec- 
tions.  The  isotherm  that  best  described  the  temperature 
characteristics  of  the  mid-waters  (40-100m)  was  chosen  from 
the  corresponding  vertical  section.   When  a  shift  was  made 
from  one  isothermal  surface  to  another,  both  surfaces  were 
drawn  to  assist  in  observing  continuity.   The  mid-waters 
were  chosen  for  two  reasons:   (1)  rapidly  changing  thermal 
conditions  in  the  surface  waters  could  be  misleading  since 
they  did  not  correlate  as  well  to  mean  density  conditions 
as  did  the  deeper  waters,  and  (2)  the  40-100m  region  should 
respond  to  offshore  currents,  upwelling,  and  wind  (to  a 
lesser  extent)  and  thus  provide  mean  isothermal  surfaces 
most  likely  to  be  related  to  changes  in  mean  current. 

2 .  Description  of  Figures  34-48 

The  comments  in  paragraph  III.  E.l  concerning  loca 
tion  of  data  points,  mean  winds,  and  mean  offshore  currents 
are  applicable  to  these  figures  as  well.   The  number  of 
temperature  values  used  to  compute  the  mean  temperatures 
for  each  block  are  listed  in  Table  I  (the  values  at  50m  are 
representative  of  those  throughout  the  mid-waters) . 
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Topographies  of  the  mean  isothermal  surfaces  were 
contoured  at  10m  intervals.   Dashed  lines  indicate  probable 
extensions  of  contours  where  data  was  not  available.   Depth 
contours  of  isothermal  surfaces  are  portrayed  only  in  regions 
of  the  Bay  where  depths  of  water  arc  greater  than  or  equal 
to  the  values  of  the  contours.   The  termination  of  a  given 
contour  is  shown  by  a  short  segment  of  the  equal  valued  bot- 
tom contour  drawn  across  it. 

The  letter  H  (high)  indicates  the  greatest  depths 
of  the  isothermal  surface  and  the  letter  L  (low)  indicates 
the  shallowest  depths. 

D.    MEAN  CURRENTS  INDICATED  BY  THE  ISOTHERMAL  SURFACES 

In  considering  the  indicated  mean  currents  it  should  be 
emphasized  that  the  results  are  based  on  the  combination  of 
data  compiled  from  1929  through  1968.   The  various  regions 
of  "highs"  and  "lows"  appearing  on  the  mean  monthly  iso- 
thermal surfaces  (Figures  34-48)  cannot  be  considered  as 
fixed.   Skogsberg  [1936]  showed  that  internal  waves  and  ad- 
vection,  acting  in  the  southern  portion  of  the  Bay,  caused 
daily  (and  even  hourly)  changes  in  the  thermal  structures. 
For  any  given  month  the  contours  vary  from  day  to  day  and 
also  from  year  to  year,  and  therefore,  the  figures  are  use- 
ful as  indicators  of  the  average  thermal  conditions  and  not 
the  conditions  for  any  given  time. 

Since  it  was  assumed  that  the  topographies  of  the  mean 
isothermal  surfaces  represented  mean  flow  (see  paragraph  C) , 
the  following  rules  apply:  (1)  mean  circulation  is  counter- 
clockwise about  regions  marked  with  an  L  (cold  water,  high 

40 


mean  density  regions) ,  and  clockwise  about  regions  marked 
with  an  H  (warm  water,  low  mean  density  regions),  (2)  mean 
flow  follows  the  contours  with  the  colder  water  always  to 
the  left,  and  (3)  stronger  currents  are  indicated  by  closer 
spacing  of  the  contours,  and  conversely. 

The  indicated  mean  currents  (inferred  from  Figures  34- 
48)  are  discussed  below.   To  prevent  repetition,  the  term 
"mean  flow"  will  be  used  to  imply  indicated  or  inferred  mean 
flow  and  the  term  "intensity"  (or  "strength")  will  imply  the 
relative  intensity  based  on  the  spacing  of  contours. 

1 .   Intensities  of  Mean  Flow  in  the  Bay  and  Offshore 
The  intensity  of  mean  flow  (I)  in  the  Bay  is  rep- 
resented by  a  formula  relating  it  to  contour  spacing  as 
follows : 

I  =  AD/5L 

where  , 

I  =  Intensity  of  mean  flow  in  mi/day, 
AD  =  The  difference  in  depth  (m)  between  contours, 
L  =  The  measured  length  (in)  between  contours. 
The  factor  of  5  (in  the  denominator)  was  applied  to  make  the 
annual  range  of  I  agree  with  the  annual  range  of  mean  off- 
shore current  (mi/day)  as  found  in  H.O.  507  (see  Figure  31). 
The  values  of  I  plotted  in  Figure  31  represent  the  maximum 
intensities  (minimum  spacing  of  contours)  for  a  given  month. 
Thus,  the  assumption  is  made  that  the  strongest  monthly  aver- 
age currents  in  the  Bay  are  equal  to  the  average  currents  off- 
shore.  Months  where  two  isothermal  surfaces  are  portrayed 
are  represented  by  an  average  value  of  I . 
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The  annual  cycle  of  the  mean  offshore  current 
(dashed  line  in  Figure  31)  was  divided  into  two  periods. 
The  California  Current  prevailed  from  February  (mid-points 
of  the  months  are  plotted  on  the  abscissa)  until  early 
November,  the  Davidson  Current  prevailing  the  remainder  of 
the  year.   The  annual  cycle  of  mean  current  intensity  in 
the  Bay  shows  peaks  in  December,  March,  and  August.   The 
peak  in  December  and  the  general  shape  of  the  curve  from 
November  to  February  is  similar  to  that  of  the  offshore 
Davidson  Current.   The  peaks  in  March  and  August  do  not  ap- 
pear to  be  related  to  peaks  of  the  offshore  current;  in 
fact,  the  shapes  of  the  Bay  and  offshore  curves  approach 
that  of  a  mirror  image  during  the  period  from  March  until 
August.   These  peaks  appear  to  h<?  related  to  upwelling  and 
are  discussed  in  subparagraph  2  below. 

2 .   Annual  Variations  of  Mean  Flow  in  the  Bay 

Inferred  mean  flow  is  depicted  by  arrows  with  the 
intensity  of  flow  (I)  indicated  to  the  nearest  whole  number 
(mi/day)  in  Figures  34-48).  ** 

a.   November  to  February  (Figures  47,  48,  and  34) 
The  variations  in  intensity  of  flow  (as  noted 
above)  are  similar  for  the  Bay  and  the  offshore  Davidson 
Current  prevailing  during  this  period.   As  described  in  para- 
graph A.l,  a  north- flowing  offshore  current  is  considered  by 
some  to  generate  a  clockwise  driven  current  within  the  Bay. 


In  this  case,  the  value  of  I  (shown  in  the  Figures) 
represent  relative  intensities  at  given  locations  and  not 
the  maximum  monthly  values . 
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Referring  to  the  mean  isothermal  surfaces,  the  current 
patterns  exhibit  this  clockwise  tendency  throughout  the 
period,  suggesting  that  the  Davidson  Current  might  be  the 
primary  generating  force. 

In  November,  the  mean  flow  enters  the  Bay  from 
the  south  and  spreads  into  the  northeast  region  with  a  broad 
clockwise  movement  indicated  in  the  center  region.   Flow 
then  exits   from  the  northern  region  of  the  Bay  (note  the 
40m  contour) . 

In  December  and  January,  the  same  general 
clockwise  tendency  of  flow  is  indicated.  However,  the  flow 
appears  to  enter  from  the  northern  region  and  exit  from  the 
southern  region  of  the  Bay.  The  difference  in  entry  points 
could  be  due  to  the  offshore  current  intensity  which  in- 
creases from  1  mi/day  in  November  to  3.4  mi/day  in  December 
and  January. 

b.   February  to  April  (Figures  35-37) 

In  February,  the  "high"  (located  near  the  cen- 
ter of  the  Bay)  intensifies  and  the  resulting  intensity  of 
flow  doubles  over  that  in  January.   Flow  appears  to  enter 
and  exit  from  the  northern  region  following  a  well  defined 
clockwise  path  throughout  the  Bay.   The  clockwise  tendency 
of  flow  could  not  have  resulted  from  offshore  current  genera- 
tion since  the  direction  of  the  offshore  current  is  away 
from  and  nearly  perpendicular  to  the  coast  (and  also  at  its 
annual  minimum  speed  of  1  mi/day) .   The  maximum  rates  of  up- 
welling  occurring  in  February  and  March  (as  discussed  in 
Chapter  III,  paragraph  B)  may  explain  the  inferred  clockwise 
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tendency  of  the  currents.   As  previously  postulated,  the 
spreading  cr  up-surge  of  the  cold  waters  from  the  Canyon 
into  the  Shallows  regions  cause  a  relatively  warm  tongue 
(high)  in  the  Canyon  region.   The  high  established  by  this 
up-surge  might  tend  to  overshadow  the  influence  of  the  weak 
flowing  California  Current  during  February. 

The  10°C  and  11°C  isothermal  surfaces  were  con- 
structed for  March.   Both  have  similar  patterns  but  the 
gradients  are  stronger  at  the  greater  depths  (10°C  surface, 
Figure  37) .   This  further  indicates  the  influence  of  up- 
welling  since  upwelling  rates  are  maximum  from  60-100m.   The 
currents  indicated  by  the  11°C  surface  (Figure  36)  enter  the 
Bay  from  the  south  and  follow  a  snake-like  pattern  before 
exiting  from  the  northern  region.   The  10°C  surface  (Figure 
37)  exhibits  the  same  general  isothermal  structure  except 
there  is  another  well  formed  low  observed  in  the  northern 
region  and  the  overall  intensity  of  flow  is  greater.   The 
development  of  the  "lows"  in  March  could  be  caused  by  an 
influence  of  the  California  Current  which  had  increased  to 
5  mi/day  (near  its  annual  maximum  intensity).   The  influence, 
by  the  theory  of  offshore  generation,  would  be  to  cause  a 
counterclockwise  gyre  within  the  Bay  which  would  agree  with 
formation  of  a  "low."   This  was  not  indicated  to  be  a  major 
feature  by  the  isothermal  topographies. 

c.   August  to  October  (Figures  42-44) 

This  period  is  discussed  next  because  the  mean 
flow  in  the  Bay  is  believed  to  be  a  result  of  conditions 
similar  to  those  described  above.   August  and  September  (the 
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months  of  Maximum  rates  of  downwelling)  showed  the  same  dip- 
ping of  isotherms  across  the  Canyon  and  the  same  relative 
warm  spots  in  the  Canyon  region  as  were  observed  during 
February  and  March.   The  possible  relationship  was  discussed 
in  Chapter  III,  paragraph  B.   The  well  established  highs  in 
August  and,  for  the  10°C  surface,  in  September  (the  deeper 
surface  where  downwelling  rates  are  also  greater)  demon- 
strate the  possible  influence  of  downwelling.   Current  trends 
are  clockwise  (as  in  February  and  March)  and  thus  opposite  to 
the  currents  theoretically  predicted  due  to  influence  by  the 
California  Current.   On  the  other  hand,  the  11°C  surface 
(shallower  waters)  in  September  exhibits  more  complex  trends 
that  could  have  been  partly  influenced  by  the  California 
Current.   Current  patterns  are  counterclockwise  in  the  mid- 
dle and  northern  regions  and  a  clockwise  trend  is  observed 
in  the  southern  region.   Downwelling  is  weaker  above  50m 
(in  September)  and  apparently  does  not  exert  as  great  an 
influence  on  the  indicated  currents.   It  was  also  observed 
that  flow  indicated  by  the  shallow  waters  (11°C  surface)  is 
nearly  opposite  in  direction  to  that  indicated  by  the  deeper 
waters  (10°C  surface)  .   Since  the  whole  region  of  the  mid- 
waters  (40-100m)  is  assumed  to  represent  the  indicated  cur- 
rents, it  is  concluded  that  the  shallow  and  the  deep  surfaces 
(in  September)  combined  in  some  fashion  to  represent  mean  flow 
The  best  approximation  of  the  mean  flow  is  thus  a  combination 
of  flows  inferred  from  both  surfaces. 
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d.  October  (Figures  45  and  46) 

The  current  trends  indicated  by  the  11°C  iso- 
thermal surface  are  similar  to  those  of  the  11°C  surface  in 
September  (Figure  44) .   Flow  enters  the  southern  region  of 
the  Bay,  follows  an  "S- shaped"  path,  and  exits  from  the 
northern  region.   The  high  observed  in  the  southern  region 
appears  to  be  a  holdover  from  the  high  observed  during  July 
and  August. 

The  12°C  surface  (Figure  46)  is  generally 
similar  to  the  11°C  surface  except  the  high  is  absent.   The 
general  counterclockwise  trends  (expecially  for  the  12°C 
surface) conform  to  the  theory  of  generation  by  an  offshore 
current . 

In  general,  it  appeared  that  weaker  offshore 
currents  came  into  the  Bay,  while  stronger  currents  caused 
a  reverse  eddy  within  the  Bay. 

e.  April  to  August  (Figures  38-41) 

The  intensity  of  indicated  mean  flow  is  low 
throughout  this  period.   Figure  31  shows  the  weak  mean  cur- 
rents which  decrease  to  an  annual  minimum  in  June. 

Current  trends  in  April  are  similar  to  those 
in  March  (Figures  36-37)  but  the  intensity  of  flow  is  only 
half  as  great.   The  weakening  of  the  high  in  the  center 
region  corresponds  to  the  decrease  in  the  rate  of  upwelling 
which  ends  altogether  in  May  to  June.   The  low  in  the  southern 
region  remained  the  same  strength  as  it  was  in  March  but  in- 
truded further  into  the  Bay  (as  noted  by  the  40m  contour  in 
Figure  38) . 
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The  indicated  mean  flow  during  May  is  almost  due 
south,  dipping  into  the  Bay  and  broadly  following  the  shore- 
line.  The  direction  of  flow  corresponds  closely  to  the  di- 
rection of  mean  offshore  current  in  the  central  and  northern 
regions  and  tends  seaward  in  the  southern  region. 

A  high  develops  in  June  which  extends  through- 
out the  Bay.   Mean  flow  enters  the  Bay  from  the  north  and 
follows  a  broad  clockwise  path,  exiting  also  from  the  north- 
ern region.   Intensity  of  mean  flow  is  at  its  annual  minimum 
in  June. 

The  high  persists  in  July  and  a  low  is  seen  to 
develop  in  the  northern  region.   Flow  enters  the  Bay  from 
the  north  and  is  observed  to  exit  from  the  southern  region. 

The  development  of  the  high  in  June  and  its 
subsequent  strengthening  from  June  through  August  (Figure 
42)  could  have  been  partly  due  to  downwelling  which  com- 
mences in  June  and  ends  in  early  November  (see  Chapter  III). 

In  summary,  the  inferred   mean  flow  in  the  Bay 
is  based  solely  on  the  topographies  of  mean  monthly  iso- 
thermal surfaces.   Inferred  current  trends  could  be  related 
to  generation  by  offshore  currents  only  during  the  Davidson 
Current  Period  and  possibly  in  the  month  of  May.   Trends  of 
inferred  current  during  the  remainder  of  the  year  appear  to 
be  influenced  more  by  up/downwelling  than  by  offshore  cur- 
rents.  The  possibility  of  up/downwelling  exerting  a  greater 
influence  than  offshore  currents  on  mean  flow  is  difficult 
to  test  and  verify.   However,  the  maximum  intensity  of  mean 
offshore  current  during  any  month  of  upwelling  or  downwelling 
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is  only  5  ni/day  (about  l/5th  of  a  kt) ,  and  the  annual  maxi- 
mum reaches  only  6.7  mi/day.   Thus,  it  would  appear  that 
the  relatively  small  horizontal  density  changes  (assumed  to 
be  mainly  a  result  of  the  surging  of  the  cold  waters  during 
up/downwelling)  between  the  waters  of  the  Canyon  and  those 
of  neighboring  shallow  areas  could  be  equally  (if  not  more) 
responsible  for  influencing  mean  Bay  currents.   The  above 
phenomena  are  of  course,  related  in  some  complex  manner. 

There  is  insufficient  data  concerning  pre- 
viously measured  currents  in  the  Bay  to  verify  the  influence 
of  up/downwelling  on  Bay  currents.   However,  the  author 
found  several  limited  current  studies  which  could  be  com- 
pared to  the  inferred  Bay  currents  of  this  report. 

E.    COMPARISONS  OF  INFERRED  MEAN  CURRENTS  TO  PREVIOUSLY 
MEASURED  CURRENTS 

Comparisons  between  the  inferred  mean  monthly  flow 
(based  on  many  years  of  data)  and  several  short  term  (6-48 
hours)  current  measurements  during  a  given  month  cannot  be 
expected  to  show  more  than  broad  agreement  for  specific 
times.   The  positions  of  the  highs  and  lows  (and  hence  the 
distribution  of  indicated  current)  on  the  mean  isothermal 
surfaces  for  a  given  month  represent  the  most  likely  con- 
ditions over  a  long  time  period  and  not  those  for  a  given 
day  of  the  month,  or  even  for  a  given  year  for  that  month. 
Finally,  the  short  term  currents  must  be  considered  only  as 
general  trends  since  their  directions  and  speeds  are  heavily 
influenced  by  tidal  forces  and  local  winds.   The  following 
comparisons  are  more  encouraging  than  the  above  limitations 
suggest.  4g 


1 .  March 

Stevenson  [1964]  deployed  three  surface  drogues  in 
early  March  from  a  location  1/2  mi  offshore  in  the  southern- 
most area  of  the  Bay.   All  drogues  moved  shoreward  (until 
beaching)  as  shown  by  the  dotted  arrows  in  Figure  36.   The 
shoreward  movements  of  the  drogues  would  appear  to  agree 
with  the  inferred  currents  which  tend  shoreward  in  the 
region  where  the  drogues  were  deployed. 

2 .  July 

McKay  [1970]  conducted  current  measurements  in  late 
July  using  a  Geomagnetic  Electrokinetograph  (G.E.K.).   The 
measured  currents  (indicated  by  dotted  arrows  in  Figure  41) 
did  not  fully  agree  with  the  plotted  indicated  current  pat- 
terns: however,  the  position  of  the  "low"  in  the  northern 
region  would  support  such  current  movement.   The  shape  of 
the  "low"  could  not  be  determined  since  there  was  no  data 
for  adjacent  blocks  to  the  north  and  east. 

3 .  August 

McKay  conducted  current  measurements  with  the  G.E.K 
in  early  and  mid-August,  1969  (solid  short  arrows  in  Fig- 
ure 49)  and  Stoddard  [1971]  employed  surface  drogues  in  mid- 
August,  1970  (dotted  arrows). 

a.   McKay's  G.E.K.  Measurements 

The  shorter  solid  arrows  represent  individual 
G.E.K.  measurements.   Considering  the  region  bounded  by  the 
60m  and  80m  contours,  the  correlation  of  the  measured  cur- 
rents to  the  inferred  mean  currents  is  very  close,  with 
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about  3/4ths  of  the  G.E.K.  measurements  observed  to  follow 
the  contours.   In  the  region  bounded  by  the  80m  contour,  the 
measured  currents  broadly  exhibit  a  clockwise  motion  with  a 
general  confusion  in  direction  noted  (which  agrees  with  the 
plotted  position  of  the  high) . 

The  G.E.K.  measurements  to  the  south  of  the  60m 
contour  appear  somewhat  confused  in  direction  when  considered 
as  a  whole.   The  arrows  with  tails  represent  measurements 
taken  on  5  August  while  all  other  solid  arrows  represent 
measurements  taken  from  11-12  August.   During  11-12  August 
the  measured  currents  in  the  extreme  southern  region  are 
seen  to  move  northward  and  then  to  turn  towards  the  west. 
This  trend  would  appear  to  agree  with  the  mean  flow  which 
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southern  region  and  could  force  the  northerly  current  to 
turn  westward.   The  currents  measured  on  5  August  tend  to- 
wards the  southwest  which  could  be  due  to  the  conditions 
such  as  displayed  by  the  isothermal  surfaces  in  July  (note 
the  30m  contour  in  Figure  41) .   Since  all  surfaces  are  as- 
sumed to  represent  the  mid-month  conditions,  the  mean  flow 
in  July  might  influence  that  in  early  August.   The  above 
discussion  is  conjectural  since  insufficient  data  is  avail- 
able to  construct  contours  in  the  region  of  interest.   The 
intention  is  to  shew  that  the  current  trends  by  McKay  in 
this  (southeastern)  region  show  plausible  agreement  with 
those  inferred  from  the  isothermal  surfaces.   The  measured 
current  velocities  were  all  significantly  higher  than  inferred 
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values  (I) .   This  is  as  expected  since  they  were  virtually 
instantaneous  measurements. 

b.   Stoddard's  Drogue  Study 

The  longer  dotted  arrows  in  Figure  49  represented 
surface  drogue  tracks.   Six  drogues  were  tracked  for  about  17 
hours  and  one  (the  southernmost  shorter  track)  was  tracked 
for  7  hours . 

The  two  tracks  in  the  region  bounded  by  the  70m 

contour  follow  definite  clockwise  gyres,  not  too  far  removed 

•   ■ 
from  the  location  of  the  high  (generally  'agreeing  with  its 


plotted  position)  .~:  The— remaining  tracks  do  not  correlate 
with  mean  flow;  in  fact  the  three  tracks  to  the  south  are 
nearly  180  degrees  out.   Data  "was  not  available  during  August 
in  Blocks  10,  11,  14,  and  21S  (see  Figure  1)  and  therefore, 
the  isothermal  surfaces  (as  shown  for  this  region)  could  be 
in  error. 

4 .   September 

Stoddard  tracked  six  surface  drogues  in  late 
September,  1970.   The  tracks  show  good  correlation  to  the 
shallower  11°C  isothermal  surface  (Figure  44)  but  generally 
are  not  well  related  to  the  10°C  surface  (Figure  43)  .   The 
following  discussion  is  based  on  comparison  to  the  11°C  sur- 
face (although  the  best  approximation  of  indicated  mean  flow 
would  lie  somewhere  in  between  the  two  surfaces) . 

The  duration  of  the  drogue  tracks  varied  from  12  to 
25  hours.   The  author  computed  the  speed  of  each  drogue  based 
on  the  distance  made  good  from  deployment  until  the  end  of 
tracking.   The  average  value  of  the  six  speeds  was  2.75  mi/day 
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which  compared  favorably  to  the  arbitrary  value  of  I  for 
September  of  3.6  mi/day  (Figure  31).   The  value  of  I  rep- 
resents the  maximum  mean  intensity  (minimum  contour  spacing) 
for  a  given  month,  and  was  made  to  agree  (in  total  range  of 
values)  with  the  intensity  of  mean  flow  for  the  offshore 
current  (see  paragraph  D.l). 

Drogues  1-3  (see  circled  numbers  adjacent  to  dotted 
arrows  in  Figure  44)  all  exhibit  clockwise  trends  which 
broadly  agree  with  the  location  of  the  high  in  the  south- 
eastern region  of  the  Bay.   Closer  correlation  is  not  pos- 
sible since  no  data  exists  in  the  block  to  the  right  of  the 
high  (shown  crossed  out  in  Figure  1) . 

Tracks  of  drogues  4  and  S  cut  directly  across  the 
40m  contour;  however,  their  tendencies  to  move  to  the  right 
when  seaward  of  Point  Pinos  agree  with  the  general  mean  cur- 
rent patterns  (note  the  north- flowing  mean  current  which 
would  force  the  currents  indicated  by  tracks  4  and  5  to  the 
right) . 

The  track  of  drogue  6  agrees  with  indicated  mean 
flow,  following  the  40m  contour  and  tending  in  a  slight 
clockwise  path.   The  speed  of. drogue  6  (computed  as  above) 
was  computed  to  be  4.1  mi/day  versus  an  arbitrary  value  of 
I  eotual  to  4  mi/day  in  the  same  region. 
5 .   October 

Stoddard  tracked  four  surface  drogues  in  early 
October,  1970.   The  tracks  show  excellent  correlation  to  the 
shallower  12°C  surface  (dotted  arrows  in  Figure  46)  and  fair 
correlation  to  the  11°C  surface  (Figure  45) .   Drogue  3  was 
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tracked  for  10  hours  and  all  other  drogues  were  tracked  for 
only  six  hours.   The  average  computed  speed  for  drogues  1, 
2,  and  4  was  eight  mi/day,  while  drogue  3  made  good  a  speed 
of  6  mi/day.   The  speeds  inferred  from  contour  spacing  are 
only  2  to  3  mi/day.   Since  the  duration  of  tracking  was 
short,  the  effects  of  the  tidal  currents  (and  local  winds) 
could  have  introduced  considerable  influence  on  the  speeds 
and  directions  of  the  drogue  tracks. 

All  drogue  tracks  broadly  followed  the  30m  contour 
(Figure  46)  and  all  except  drogue  3  moved  in  a  direction 
related  to  the  indicated  mean  current.   The  path  of  drogue 
4  almost  exactly  agrees  with  the  indicated  flow. 
6 .   November 

Stoddard  tracked  eight  surface  drogues  in  mid- 
November,  1970  (shown  by  solid,  dashed,  and  dotted  arrows, 
numbered  1-8  in  Figure  50).   The  duration  of  tracking  varied 
from  6-42  hours.   The  computed  speeds  (again  based  on  dis- 
tance made  good  per  total  tracking  time)  are  listed  below: 

Duration  of  Tracking    Computed  Speed 
Drogue  No.  (hr)  (mi/day) 

1  .42  5.7 

2  22  3.9 

3  21  1.2 

4  23  3.4 

5  7.5  11.2 

6  6.5  14.1 

7  6.5  '  10.7 

8  6  10.8 

The  consistently  greater  speeds  for  drogues  5-S 
were  probably  due  to  tidal  currents  prevalent  in  the  short 
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duration  of  tracking  time  and  were  thus  disregarded.   Drogues 
1-4  exhibited  speeds  closely  related  to  the  arbitrary  value 
of  I  for  November.   The  combined  average  speed  of  drogues 
1-4  was  3.6  mi/day  while  the  arbitrary  maximum  value  of  I 
was  4.2  (Figure  31).   The  average  value  of  I  for  November 
would  be  somewhat  lower  (but  nearly  impossible  to  estimate 
from  contour  spacing  alone) . 

The  tracks  of  all  drogues  exhibited  clockwise 
tendencies  which  correlated  with  indicated  mean  current 
about  the  extensive  high  in  the  center  of  the  Bay.   The  in- 
dividual tracks  of  drogues  1  and  8  closely  approximate  the 
indicated  mean  current  while  the  remaining  drogues  showed 
fair  correlation  (Figure  50) . 

It  is  noted  that  the  speed  of  drogue  3  (1.2  mi/day) 
is  significantly  lower  than  the  speeds  of  drogues  1,  2,  and 
4.   Also,  the  track  of  drogue  3  describes  a  tight  clockwise 
path  which  is  not  the  case  for  the  other  tracks.   This  led 
the  author  to  believe  that  drogue  3  might  have  been  nearer 
the  center  of  a  high  than  is  shown  in  Figure  50.   It  is  as- 
sumed that  the  high  (and  the  entire  mean  isothermal  surface 
as  well)  could  probably  vary  in  location  from  day  to  day 
(and  year  to  year)  for  a  given  month.   On  this  basis,  the 
mean  isothermal  surface  for  November  was  translated  to  a 
position  where  the  high  overlayed  the  track  of  drogue  3. 
The  distance  of  translation  was  3  mi  in  a  west- southwest 
direction.   Figure  51  shows  the  results  of  the  translation. 
All  drogue  tracks  are  now  observed  to  closely  follow  the 
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translated  contours  and  the  correlation  to  indicated  mean 
current  is  excellent. 

F .    SUMMARY 

The  general  agreement  between  measured  currents  from 
past  studies  and  indicated  mean  currents  for  this  study  is 
encouraging.   The  total  number  of  comparisons  made  is  far 
from  that  required  to  fully  substantiate  the  inferred  mean 
currents  of  this  study;  however,  the  comparisons  do  lend 
considerable  support  to  the  feasibility  of  using  mean  iso- 
thermal surfaces  as  indicators  of  mean  flow,  which  was  a 
primary  objective  of  this  study. 


V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

1.  Progressive  warming  with  time  is  indicated  through- 
out Monterey  Bay,  mainly  occurring  at  depths  from  the  surface 
to  300m.   The  temperature  was  noted  to  increase  as  much  as 
1.8°C  over  the  29  year  interval  from  1931  to  1960. 

2.  The  intensity  of  upwelling  varies  in  different  geo- 
graphic regions  of  the  Bay.   This  variation  in  the  maximum 
rates  of  upwelling  (and  downwelling)  is  believed  to  be  a 
major  factor  controlling  the  shapes  of  mean  isothermal  sur- 
faces (and  indicated  mean  currents  as  well)  during  the  months 
of  maximum  up/ downwelling . 
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indicators  of  mean  current  shows  promise.   The  comparisons 
of  inferred  mean  flow  to  direct  current  measurements  (as 
shown  in  Figures  34-48)  demonstrates  general  agreement  be- 
tween the  two. 

B.  RECOMMENDATIONS 

1.   The  data  collected  throughout  Monterey  Bay  is  in- 
sufficient for  several  of  the  blocks  in  the  author's  grid 
system  (Figure  1).   To  obtain  a  more  complete  picture,  both 
in  the  mean  monthly  sea  surface  temperature  distribution 
and  in  the  topographies  of  mean  monthly  isothermal  surfaces , 
a  continued  search  of  data  for  the  following  blocks  would  be 
beneficial:   7,  8,  9,  10,  11,  14,  16,  18,  21S  as  well  as  the 
block  shown  crossed  out  in  Figure  1. 
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2.  Future  studies  should  be  made  of  the  mean  monthly- 
current  distribution  in  the  Bay  (the  author  was  unable  to 
locate  any  previous  work  in  this  area) .   Whatever  methods 
are  used  in  the  future  might  be  compared  to  the  mean  flow 
inferred  in  this  study  for  purposes  of  attempted  verifica- 
tion. 

3.  To  aid  in  the  future  study  of  oceanographic  para- 
meters (in  general)  of  Monterey  Bay.,  a  computer  storage  and 
retrieval  system  should  be  instituted  and  maintained  on  a 
continuing  basis.   An  example  of  a  computer  system  designed 
for  this  purpose  is  presented  in  Appendix  B  (the  "final 
computer  program,  number  11"). 

Closer  coordination  with  NODC  is  required  of  all 
oceanographic  facilities  so  that  the  oceanographic  data 
collected  in  Monterey  Bay  over  the  years  can  be  made  avail- 
able to  any  investigator  desiring  its  use. 
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APPENDIX  A:      DATA   SOURCES 

A.  PROBLEMS  ENCOUNTERED 

The  only  previous  research  found  that  categorized  oceano- 
graphic  data  collected  in  Monterey  Bay  was  that  presented  in 
an  NPS  research  paper  by  P.  M.  Magruder,  Jr.,  [1969]. 
Magruder 's  paper  was  most  beneficial  in  delineating  the  major 
sources  of  data.   Many  reports  were  located  covering  related 
subjects   such  as  biology,  sound  velocity  variation,  and 
light  attenuation  studies.   Many  of  these  reports  used  tem- 
perature and  salinity  data.   However,  the  data  was  generally 
not  contained  therein,  nor  is  there  any  reference  to  its  lo- 
cation in  some  reports.   There  is  a  general  lack  of  continu- 
ity in  reporls  which  contained  temperature  and  salinity  ^?ta, 
Methods  of  collection  were  sometimes  omitted  and  very  seldom. 
was  there  any  information  describing  environmental  condi- 
tions, etc. 

B.  MAJOR  DATA  SOURCES  (FIGURE  3) 
1.   Bolin  Data 

Rolf  L.  Bolin  and  collaborators  [1964]  collected 
Nansen  samples  over  a  five  year  period  from  1951-1955  at  the 
position  36-42N,  122-02W.   Casts  of  24-bottles  were  gener- 
ally employed  to  depths  of  1000m  at  weekly  intervals.   A 
total  of  233  casts  were  made  and  analyzed  for  temperature, 
salinity,  and  other  chemical  parameters.   A  copy  of  this  re- 
port is  on  file  at  the  Naval  Postgraduate  School  Library. 
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2 .  Fleet  Numerical  Weather  Center  (FNWC)  XBT  Data 
Nansen  cast,  bathythermograph  (BT)  ,  experimental 

BT  (XBT),  and  message  BT  data  are  on  file  at  FNWC.   The 
Nansen  cast  and  BT  data  are  also  filed  at  the  National 
Oceanographic  Data  Center  (NODC)  and  are  kept  more  com- 
pletely there. 

FNWC  is  presently  digitizing  the  XBT  data  for  in- 
clusion in  computer  tape  files  at  NODC.   Two  9 -track  com- 
puter tapes  (NPS  264  and  NPS  265),  held  by  Dr.  Denner  at 
the  Naval  Postgraduate  School  Oceanography  Department, 
contain  the  world  oceans  XBT  data  which  consists  of  about 
40,000  individual  stations.   A  computer  program  (see  article  B.2, 
Appendix  B)  written  by  the  author  reveals  only  4  3  XBT  sta- 
tions scattered  randomly  throughout  the  Bay.   There  have 
been  considerably  more  than  43  XBT  drops  in  the  Bay  over  the 
last  three  years  but  the  data  is  not  as  yet  available. 

3.  Hcpkins-CalCOFI  Data 

Hopkins  Marine  Station  has  conducted  a  continuous 
hydrobiological  survey  of  Monterey  Bay  for  the  California 
Cooperative  Oceanic  Fisheries  Investigations  (CalCOFI)  since 
1951.   Commencing  in  March  of.  1954  the  data  (unpublished) 
was  tabulated  and  distributed  in  annual  reports  containing 
monthly  averages  of  Bucket,  Nansen,  and  BT  temperatures, 
salinities,  and  densities  (sigma-T)  for  six  oceanographic 
stations  in  the  Bay.   The  averages  were  based  on  approxi- 
mately weekly  visits  to  each  station;   [Hopkins  1954-1967]. 
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a.  Ocean  Temperatures 

Surface  temperatures  were  read  from  a  bucket 
thermometer  and  BT  data  was  listed  at  10m  intervals  from 
the  surface  to  50m  for  stations  2,  3,  4,  and  6,  and  to  30 
and  20m  for  stations  1  and  5,  respectively.   This  data  was 
not  computerized  since  the  majority  of  it  was  included  as 
NODC-BT  data  in  paragraph  4  below. 

b.  Salinity  and  Sigma-T  Data 

Nansen  samples  taken  at  the  surface  and  15m 
were  titrated  for  chlorinity  and  thence  salinity  and  sigma-T 
were  determined.   Temperature,  salirity,  and  sigma-T  were 
tabulated  along  with  their  monthly  averages  for  each  station 
4 .   National  Oceanographic  Data  Center  (NODC)  Data 
Magruder  obtained  complete  printed  listings  of 
Nansen  cast  and  BT  data  from  NODC  ir.  January  of  1969.   The 
area  of  cov3rage  included  Marsden  Squares  121  (61  and  62). 
The  listings  are  available  in  the  NPGS  Oceanography  files. 

For  purposes  of,  (1)  rapid  computer  processing, 
and  (2)  updating  of  data  through  1970;  magnetic  tapes  were 
requested  from  NODC.   Two  7- track  tapes  were  received  in 
December  1970,  one  listing  Nansen  cast  data  in  Marsden 
Square  121(R1019)  and  the  other  listing  BT  data  (R1030) . 
The  formats  of  these  tapes  were  different  and  had  to  be  pro- 
cessed separately  (see  articles  B.4,  5  of  Appendix  B) .   The 
tapes  are  the  property  of  the  NPGS  Oceanography  Department 
(see  Dr .  Leipper) . 
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The  total  number  of  Nansen  and  BT  stations  within 
Monterey  Bay  (retrieved  from  above  rapes)  were  107  and  1377 
respectively . 

a.  Nansen  Cast  Data 

The  data  included  temperature,  salinity,  oxygen, 
phosphates,  silicates,  nitrates,  and  pH  as  observed  para- 
meters.  All  values  were  interpolated  to  standard  depths  by 
a  three  point  Lagrange  method.   Each  station  is  identified 
by  a  Master  Card,  which  lists  such  information  as  ship  name, 
position,  date,  etc. 

b.  BT  Data 

BT  temperatures  were  listed  at  5m  intervals 
from  the  surface  to  maximum  depth.   Each  station  was  again 
identified  by  a  Master  Card. 

5 .  Naval  Postgraduate  School  fNPGS)  Data 

Four  stations  were  established  by  NPGS  in  October 
1966.   Nansen  casts  and  BT  drops  were  made  regularly  and 
all  chemical  parameters  were  filed  quarterly. 

A  continuing  program  of  computer  processing  and 
filing  of  data  cards  for  these  stations  is  in  progress. 

6 .  Skogsberg  Data 

Tage  Skogsberg  [1936]  conducted  an  exhaustive  study 
of  thermal  conditions  in  the  southern  portion  of  Monterey 
Bay,  collecting  over  15,000  samples  from  1929-1933.   Data 
included  temperature,  chlorinity,  phosphate,  and  silicate 
from  23  separate  stations.   The  original  data  is  on  file  at 
Hopkins  Marine  Station.   Skogsberg  tabulated  monthly  temperature 
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averages  for  two  of  his  stations ,  both  of  which  were  used  in 
the  author's  analysis  (Chapters  III  and  IV). 
7 .   Sea  Surface  Temperatures 

a.  Hopkins  Marine  Station 

Sea  surface  temperature  and  salinity  has  been 
measured  daily  at  the  station  since  1919.   All  records  are 
filed  by  Hopkins  and  monthly  averages  are  published  by 
Scripps  Institution  of  Oceanography  under  the  title: 
"Data  Report,  Surface  Water  Temperatures  at  Shore  Stations,- 
United  States  West  Coast  and  Baja  California." 

b.  Santa  Cruz 

Sea  surface  temperature  has  been  measured  daily 
since  1956  and  data  is  also  published  by  Scripps  under  the 
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C.    MINOR  DATA  SOURCES 

The  below  listed  sources  comprise  data  covering  either 
small  areas  or  that  limited  to  only  several  observations. 
A  complete  listing  was  impossible. 

1 .  "Monterey  Bay  Bibliography"  [Baron.  1971] 

This  comprehensive  study  was  recently  completed  by 
the  Moss  Landing  Marine  Laboratories'  of  the  California  State 
Colleges.   It  includes  all  known  scientific  literature  con- 
cerning Monterey  Bay.   Limited  copies  are  now  being  published 

2 .  Yeske  and  Waer  Data 

Yeske  and  Waer  [1968]  studied  light  attenuation  at 
two  stations  in  the  Bay  during  the  period  26  July  to  31 
August,  1968.   A  total  of  39  Nansen  casts  were  made  and  BT 
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drops  accompanied  each  cast.   Data  was  filed  with  NODC  and 
is  contained  on  the  computer  tapes  referenced  in  para- 
graph C.4  above. 

3.   Message  BT  Data 

Very  little  data  was  found  at  FNWC  (only  21  sta- 
tions for  all  of  Marsden  Square  121)  and  none  more  recent 
than  1968. 

4 ■   OC-3520  Data 

The  Chemical  Oceanography  Course  (0C-   3520)  at  NPGS 
sponsors  limited  duration  Bay  cruises  for  students.   Various 
oceanographic  parameters  are  collected  and  final  reports 
are  required.   The  data  is  not  currently  being  filed  but 
some  of  it  is  available. 

D.    DISCUSSION 

It  was  clear  that  a  complete  source  listing  of  the 
chemical  parameters  in  the  Bay  would  be  a  monumental  task. 
NODC  maintains  a  computer  listing  of  oceanographic  data 
throughout  the  world's  oceans  but  little  of  the  data  col- 
lected in  Monterey  Bay  has  ever  reached  there.   Only  107 
Nansen  casts  were  recorded  by  NODC  for  the  last  40  odd  years, 
representing  a  trifle  of  the  data  actually  collected. 
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APPENDIX  B:   COMPUTER  PROCESSING  OF  DATA 

Two  major  goals  were  sought  in  designing  the  "final 
computer  program":   (1)  to  store  temperature,  salinity,  and 
density  data  in  a  form  suitable  for  the  computation  of  mean 
monthly  averages  at  given  geographic  locations,  and  (2)  to 
design  a  flexible  system  for  storage  and  retrieval  of  various 
oceanographic  parameters  on  a  continuing  bases. 

The  second  goal  was  achieved  by  designing  the  "station 
title  card"  in  such  a  manner  that  a  simple  extension  could 
be  made  to  the  "station  data  card"  to  include  additional 
selected  oceanographic  parameters  without  modification  of 
existing  formats  (see  paragraphs  A. 2. a  §  b) .   The  institution 
ui  sUl-Ti  a  caiva  DV.5L.sru  cg-uj.u  maicc  one  re  tncvai  an*j.  usage  Cj. 
oceanographic  data  from  Monterey  Bay  a  very  simple  task. 

A.   STORAGE  OF  TEMPERATURE,  SALINITY"  AND  DENSITY  (SIGMA-T) 
DATA 

1 .   Problems  Encountered 

About  60%  of  available  data  was  already  in  computer 
format,  either  on  magnetic  tape  or  on  data  cards.   Regret- 
tably, no  two  formats  were  consistent.   Thus,  a  separate 
computer  program  had  to  be  written  to  reduce  each  data  source 
to  the  author's  format.   This  often  included  calculation  of 
densities  (when  they  were  not  included  as  part  of  the 
original  format)  and  always  required  interpolation  of  one 
form  or  another  to  place  the  data  at  the  author's  selected 
"standard  depths"  (standard  depths  are  listed  in  Table  IV) . 
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The  reduction  of  data  already  i/i  conputer  format  was  actu- 
ally more  time  consuming  than  processing  the  tabulated  data. 

The  remainder  of  the  available  data  was  tabulated 
in  the  various  reports  cited  in  Appendix  A.   First  it  was 
keypunched  onto  computer  cards  in  various  formats,  depending 
on  the  format  of  the  data  source.   The  next  step  was  the 
printing  out  of  data  on  computer  paper  to  compare  it  for 
accuracy  against  the  original  source.   After  proofreading, 
computer  programs  were  written  to  convert  the  data  on  punched 
cards  to  the  desired  format.   Again,  this  included  calcula- 
tions and  interpolations  for  the  desired  results. 

The  last  step  in  the  process  was  to  check  the 
various  outputs  for  correctness  and  then  write  the  data  onto 
a  computer  "master  tape."   Paragraph  B  describes  the  various 
computer  programs  to  the  extent  of  including  information  not 
apparent  from  inspection  of  the  programs  themselves. 
2 .   Format  of  the  Master  Tape 

The  master  tape  (NPS  serial  number  306)  served  as 
the  storage  location  for  all  data  processed.   The  general 
format  of  data  includes  a  station  title  record  (identifying 
the  oceanographic  station)  and  a  station  data  record  (list- 
ing the  oceanographic  parameters).   Table  II  lists  all  data 
contained  on  the  master  tape. 

a.   The  Station  Title  Record 

The  station  title  record  allows  the  flexibility 
required  in  general  retrieval  of  data.   Referring  to  Table 
III,  the  "80-column  record  numbers"  refer  to  the  inclusive 
positions  on  the  station  title  record  where  the  information 
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TABLE  II.   Listing  of  Specific  Data  Contained  on  Computer 
Master  Tape,  NPS  506 


Data  Source Data Year  (incl .) Remarks 

Bolin  T,S,at       1951-1955     229  Nansen  stations 

FNWC  T        Various      40  XBT  stations 

Hopkins-CalCOFI   T,S,at       1954-1967     972  Nansen  stations, 

each  station  repre- 
senting the  monthly 
average  of  4  stations 
at  0  and  15m. 

NODC  See  Chap  II  Various      107  Nansen  stations 

NODC  T        Various      1377  BT  stations 

NPGS  See  Chap  II  Sept  1966-    175  Nansen  stations 

Sept  1967 


Note  1:   See  Appendix  A  for  detailed  discussion  of  data 
sources . 

Note  2:   In  cases  where  only  portions  of  a  data  source  were 
computerized,  the  year  and  type  of  data  defines 
such  portions. 

Note  3:   The  "Data  Source"  is  listed  by  author  or  collecting 
activity . 
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TABLE  III.   The  Station  Title  Record 


80-Column 

Record  No.    Variable 


Definition 


1-2 

NCARDS 

4-7 

NANS,  XBT, 
BT,  STD 

9-10 

NDAY 

12-13 

NMO 

15-16 

NYEAR 

18-21 


43-45 


48-63 


NTIME 


24- 

-25 

LATDEG 

27- 

-30 

RLATMI 

33- 

-35 

LONDEG 

37- 

-40 

RLONMI 

NWF 


The  total  number  of  station  data 
records  following  the  title  record 

A  left-adjusted  variable  defining 
the  type  of  data  collected. 

The  day  data  collected. 

The  month  data  collected. 

The  last  two  digits  of  the  year 
data  collected. 

The  time  (GMT)  data  collected.   If 
unknown,  NTIME  appears  as  0000. 

The  degrees  of  station  latitude. 

The  minutes  and  tenths  of  minutes 
of  station  latitude. 

The  degrees  of  station  longitude. 

The  minutes  and  tenths  of  minutes 
of  station  longitude. 

NWF=999  specifies  Hopkins -CalCOFI 
monthly  average  values  (4  stations 
per  average) . 

Ship  name. 


or 


48-49 
50-51 
52-55 
56-59 


NODC  country  code. 

NODC  ship  code. 

NODC  station  I.D.  number. 

NODC  Stano  number. 


65-80 


REF 


Reference  data  by  author  or  facility 
controlling  data  collection. 
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was  entered.   "Variable"  was  the  computer  name  (applied  by 
author)  used  in  programming.   Desired  data  may  be  retrieved 
from  the  master  tape  by  any  one,  or  any  combination  of  the 
listed  variables.   For  example,  assume  an  investigator  is 
interested  in  conducting  a  study  of  experimental  (XBT)  versus 
mechanical  (BT)  bathythermograph  data  in  the  Submarine  Can- 
yon, collected  over  the  last  five  years.   By  a  series  of 
"computer  logical  IF  statements"  [see  Blatt,  1968]  all  sta- 
tion title  records  are  eliminated  except  those  that  are 
either  XBT  or  BT ,  have  geographic  positions  in  the  Canyon, 
and  fall  within  the  time  frame  (NYEAR)  of  1966  or  later, 
b.   The  Station  Data  Record 

Each  station  data  record  contains  the  standard 
depth  of  observation  (m)  plus  one  or  more  of  the  parameters 
listed  in  Table  IV  above  the  dashed  line.   The  oceanographic 
parameters  below  the  dashed  line  were  not  included  on  the 
master  tape  but  are  listed  as  the  recommended  additions  re- 
ferred to  in  the  second  paragraph  of  this  Appendix. 

B.    COMPUTER  PROGRAMS  USED  IN  DATA  STORAGE 

Ten  computer  programs  were  required  to  compile  the  mas- 
ter tape  and  are  briefly  described  below.   The  programs  will 
be  found  in  the  section  entitled,  "Computer  Programs." 
1.   Bolin  Data 

a.   Program  Number  One 

Cards  were  punched  from  the  tabulated  data  of 
depth,  temperature,  and  salinity.  Program  number  one  cal- 
culated densities  and  linearly  interpolated  all  data  to  the 
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TABLE  IV.   The  Station  Data  Record 


80-Column 
Record  No 


Variable 


Definition 


1-5 


7-11 
13-17 
19-25 


STDEP 


TEM 


SAL 


SGT 


The  author's  selected  standard  depths 
(m)  ,  as :  0000  .   Standard  depths  used 
were:  0,  10,  15,  20,  30,  50,  75,  100, 
125,  150,  200,  250,  300,  400,  ...  , 
1400,  1500,  1750,  2000,  2500,  3000, 
4000,  .  .  ,  6000. 

The  temperature  to  hundredths  of 
degrees  C,  as:  10.25. 

The  salinity  to  hundredths  of  parts 
per  thousand,  as:   33.65. 

Sigma-T  (density)  to  hundredths  of 
units  ,  as :  26.77. 


25-27 

XOXG 

29-31 

XP04 

33-35 

XN02 

37-39 

XN03 

41-43 

SI04 

45-47 

XPH 

49-80 



Oxygen  to  tenths  of  ml  per  1. 

P04-P  to  tenths  of  g-at  per  1. 

N02-N  to  tenths  of  g-at  per  1. 

N03-N  to  tenths  of  g-at  per  1. 

Si04-Si  to  tenths  of  g-at  per  1. 

pH  of  seawater  to  tenths  of  units. 

Can  be  used  for  additional  storage 
as  required;  e.g.,  computed  sound 
velocity,  dynamic  parameters,  etc. 


Note  1:   Data  above  dashed  line  is  included  on  master  tape. 
Data  below  dashed  line  is  the  recommended  additions 
for  a  future  storage  and  retrieval  system. 
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selected  standard  depths.   Linear  interpolation  was  considered 
accurate  since  the  maximum  difference  between  tabulated  and 
standard  depths  was  only  25m.   The  final  output  consisted  of 
computer  cards  in  the  format  for  the  master  tape. 

2.  Fleet  Numerical  Weather  Center  (FNWC)  Data 
a.   Program  Number  Two 

The  existing  9-track  tapes  (NFS  264  and  265) 
containing  world  ocean's  XBT  data  through  1968  were  pro- 
cessed.  Data  in  Monterey  Bay  was  retrieved,  linearly  in- 
terpolated, and  punched  out  on  computer  cards  in  the  format 
for  the  master  tape.   Linear  interpolation  was  exact  be- 
cause XBT  data  occurred  at  every  depth  where  temperatures 
deviated  from  a  linear  path. 

3 .  Hopkins-CalCOFI  Data 

a.   Program  Number  Three 

Computer  cards  were  punched  from  the  tabulated 
annual  reports  [Hopkins  1954-1967].   Program  number  three 
processed  the  computer  cards,  linearly  interpolated  for  the 
10m  values  (a  selected  standard  depth)  ,  and  punched  new  com- 
puter cards  in  the  format  for  the  master  tape. 

Since  each  station  represented  a  monthly  average 
value  (generally  the  average  of  four  stations)  the  data  had 
to  be  processed  differently  in  the  final  program.   A  weighting 
factor  (NWF=999)  was  added  to  each  station  title  record  and 
served  as  a  branch  instruction  to  weight  this  data  "four  to 
one"  over  other  data  on  the  tape. 
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4 .  National  Oceanographic  Data  Center  (NQDC)  Nansen 
Cast  Data 

a.  Program  Number  Four 

This  program  read  data  from  a  7-track  tape 
(R  1019)  ,  computed  the  number  of  station  data  records  fol- 
lowing each  station  title  record  (NCARDS) ,  and  wrote  the 
data  (in  modified  format)  onto  a  9-track  tape  (NFS  297). 

b.  Program  Number  Five 

Data  within  Monterey  Bay  was  retrieved  from  the 
9-track  tape  (NPS  297)  and  written  directly  onto  the  master 
tape  in  the  author's  format.   One  linear  interpolation  for 
15m  was  required  per  station,  all  other  data  being  at  the 
author's  standard  depths. 

5 .  National  Oceanographic  Data  Center  (NQDC)  Bathy- 
thermograph Data 

a.  Program  Number  Six 

Data  was  read  from  a  7-track  tape  (R  10  30)  and 
written  directly  onto  a  9-track  tape  (NPS  296).   Each  sta- 
tion title  record  contained  the  number  of  following  station 
data  records,  hence,  no  further  programming  was  necessary. 

b.  Program  Number  Seven 

Data  within  Monterey  Bay  was  retrieved  from 
the  9-track  tape  (NPS  296)  and  written  directly  onto  the 
master  tape  in  the  author's  format.   No  interpolation  was 
required  as  data  was  entered  on  the  tape  (R  1030)  at  5m 
intervals . 
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6 .  Naval  Postgraduate  School  Data 
a.   Program  Number  Eight 

This  program  was  written  by  Robert  Middelberg 
for  the  Oceanography  Department.   He  was  kind  enough   to 
allow  the  author  to  copy  and  then  modify  the  program  for 
usage  herein.   Data  cards  prepared  by  Mr.  Middleberg  were 
read  by  the  modified  program  and  densities  were  computed. 
All  data  was  interpolated  to  the  author's  standard  depths 
using  Middleberg' s  "three  point  Lagrange"  method.   The 
Lagrange  interpolation  method  more  closely  approximates 
the  normal  shape  of  a  temperature  profile  when  data  is 
widespread  (50-150m  between  observations)  as  was  the  case 
here.   Finally,  data  was  written  directly  onto  the  master 

Uclpt      XJi      Lilt      ctUuilui       o      iUiiliai.. 

7 .  Transferring  Punched  Card  Data 
a.   Program  Number  Nine 

Program  number  nine  read  the  punched  output 
from  programs  1,  2,  and  3  above  and  wrote  it  onto  the  mas- 
ter tape. 

8 .  Master  Tape  Record  Count 
a.   Program  Number  Ten 

When  working  with  computer  tape  it  is  essential 
to  keep  an  accurate  record  count  of  data  as  it  is  stored; 
otherwise,  errors  that  may  be  made  would  be  nearly  impos- 
sible to  locate.   Program  number  ten  was  used  to  keep  a 
running  account  of  data  stored  and  a  double  check  for  cor- 
rect format. 


72 


C.    THE  FINAL  COMPUTER  PROGRAM 

Program  number  11  computes  the  mean  monthly  temperatures 
and  densities  (at  author's  standard  depths)  throughout  the 
Bay,  using  data  from  the  master  tape  as  its  input.   Table 
V  defines  the  output  as  it  appeared  on  the  master  printout. 
Figure  2  shows  the  initial  field  of  blocks  used  for  the 
segregation  of  data.   The  numbered  blocks  on  the  Figure  were 
the  only  ones  containing  significant  data  (at  least  one 
month  with  three  or  more  temperatures  in  the  block) .   This 
program  will  also  be  found   in  the  section  entitled,  "Com- 
puter Programs." 
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TABLE  V:   Definition  of  Computer  Variables  in  the  Final 
Output  of  Program  Number  11 


Variable 


Definition 


TEMP (AV) 

SIGMA-T(AVG) 
NR  WTD  VALUES 
SUM  WTD  TEMPS 
STDEV  WTD  TEMPS 

SUM  WTD  SIGMA-T 
STDEV  WTD  SIGMA-T 
NR  UNWTD  TEMPS 

SUM  UNWTD  TEMPS 


The  average  of  all  temperatures  (weighted 
and  unweighted)  for  a  given  block,  month, 
and  depth. 

As  above  for  density  (sigma-T) . 

The  number  of  weighted  stations  processed 

The  sum  of  all  weighted  temperature  value: 

The  computed  standard  deviation  of  all 
weighted  temperature  values. 

As  above  for  density. 

As  above  for  density. 

The  number  of  non-weighted  temperature 
values . 

Th^  cutti  nf  ?11  p on - vre i "b "*" e d  temperature 
values . 


STDEV  UNWTD  TEMPS  The  computed  standard  deviation  of  all 

non-weighted  temperature  values. 


NR  UNWTD  SIGMA-T 
SUM  UNWTD  SIGMA-T 
STDEV  UW  SIGMA-T 
WTD  TEMPS 

WTD  SIGMA-T 
UNWTD  TEMPS 

UNWTD  SIGMA-T 


As  above  for  density. 

As  above  for  density. 

As  above  for  density. 

A  listing  of  all  weighted  temperature 
values . 

A  listing  of  all  weighted  sigma-T  values. 

A  listing  of  all  non-weighted  temperature 
values . 

A  listing  of  all  non-weighted  sigma-T 
values . 


Note  1:   All  variables  above  are  for  a  given  block,  month, 
and  standard  depth. 
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NOTE:   Only  the  numbered  blocks  contained  significant  data 


figure  2.   The  Initial  Computer  Grid  Systei 
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Figure  3.   Location  of 
Major  Data  Sources. 
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TIME  ( MONTHS) 
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1000- 


P'igure  k.      Average  annual  temperature  (*C) 
cycle  based  on  five  years  of  data  (unT>er  100m 
from  1929-1933 t  below  100m  from  1932-1933)  by 
Skogsberg.   Dashed  lines  note  scale  change. 
Note:  The  9° 0  isotherms  by  Skogsberg  (S) , 
Bolin  (B)  and  CalCOFI  surveys  (C)  demonstrate 
the  downward  migration  of  isotherms. 
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NOTE:   Dashed  line  denotes  scale  change. 


Figure  5. 


Annual  .'lean  Temperature  Cycle  L:  Jolin 
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Figure  5 


Annual  Mean  Temperature  Cycle  Based  on  Ca 
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Figure  6A.   The  Progressive  '.varming  of  Bay  Waters 
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Figure  7.   Standard  Deviation  of  Temperature 
for  Selected  Blocks  as  a  Function  of  Time. 
Sea  Surface  Temperature. 
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Figure    8.      Standard  Deviation  of  Temperature, 
Averaged   in  Typical  Blocks,    as   a  Function  of 
Depth. 
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Figure    9.      The   Annual     lean  Tenperauure    L) 
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Figure  10 


The  Annual  Mean  Temperature  Cycle  at  Block  3 
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Figure  13. 


The  Annual  Cycles  of  Mean  Sea  Surface  To       ire 
in  Very  Shallow  Regions  of  the  Bay. 
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Figure  14.   The  Annual  'lean  Temperature  Cycles  at  Blocks 
1  and  5 . 
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Figure  15. 


The  Annual  Cycles  of  'lean  Sea  Surface  Te       re 
in  Relatively  Deeper  Legions  of  the  Sh 
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NOTE  : 


The  steep  slopes  of  the  10°C  and  11°C  isotherms  at 
Block  4  illustrate   the  up -slope  entry  or  cola 


waters  into  the  region 


Figure  16. 


A  s  Annual  Mean  Temperature  Cycles  at 
2  and  I  . 
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^i-yre  1"   Typical  Annual  Mean  Temperature  Cycles  of 
Canvon  versus  an  Offshore  Station. 
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Figure  19.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

JANUARY 
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Figure  20.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

FEBRUARY 
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Figure  22.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

APRIL 
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Figure  24.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

JUNE 
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Figure  25.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

JULY 
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Figure  26.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

AUGUST 
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Figure  28.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

OCTOBER 
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Figure  29.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

NOVEMBER 
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Figure  30.   Distribution 
of  the  Mean  Monthly  Sea 
Surface  Temperature  for 

DECEMBER 
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dashed  line  denotes  the  mean  offshore  flow. 
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NOTE:   The  dashed  lines  denote  mid-values  of  mean  density. 
Figure  32.   The  Annual  ilean  Density  Cycles  at  Selected  Blocks 


10  7 


TIME    (MONTHS) 


M        A 


*Mzn*^*dL  .«3Bcn«rT3gny^-.jeui  m  imli  iMW iMpww  i  i  JUnKWB 


J    J    A 


0   N   D 


DEPTH 
(m) 


1000 


500 


1000 


NOTt :   The  non-labeled  dashed  lines  denote  mid-values  of 
mean  density. 


Figure  33.   The  Annual  Mean  Temperature  versus  Dens: 


Cycles  at  Block  13 
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Figure  34.   The  Topography 
of  the  11°C  Mean  Isother- 
mal Surface  for  January. 
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Figure  35.   The  Topography 
of  the  11°C  Mean  Isother- 
mal Surface  for  February. 
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Figure  36.   The  Topography 
of  the  11 °C  Mean  Isother- 
mal Surface  for  March. 
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Figure  37.   The  Topography 
of  the  10°C  Mean  Isother- 
mal Surface  for  March. 
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Figure  40.   The  Topography 
of  the  10°C  Mean  Isother- 
mal Surface  for  June. 
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Figure  41.   The  Topography 
of  the  10°C  Mean  Isother- 
mal Surface  for  July. 
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Figure  42.   The  Topography 
of  the  10°C  Mean  Isother- 
mal Surface  for  August. 
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Figure  43.   The  Topography 
of  the  10°C  Mean  Isother- 
mal Surface  for  September. 
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Figure  44.   The  Topography 
of  the  11°C  Mean  Isother- 
mal Surface  for  September. 
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Figure  47.   The  Topography 
of  the  12°C  Mean  Isother- 
mal Surface  for  November. 
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Figure  48.   The  Topography 
of  the  12°C  Mean  Isothermal 
Surface  for  December. 
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Figure    50.    Comparison   of    In- 
ferred  Currents    to   Past   Mea- 
sured  Currents    in  November 
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